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ABSTRACT 

A  parametric  study  of  pressurized  light  water  reactor 
cores  was  undertaken,  in  which  "both  heat  transfer  and  nuclear 
influences  of  the  basic  parameters  on  the  quantity  and  required 
enrichment  of  the  uranium  fuel  were  to  be  considered. 

Methods  were  developed  for  the  following; 

1.  Computation  of  the  minimum  core  size  for  a  specified 
power  output  on  the  basis  of  heat  transfer  characteristics  of 
the  unit  cell  parameters. 

2.  Computation  of  the  critical  mass  of  U-235  required 
for  a  given  reactor  by  an  adaptation  of  a  digital  computer 
program  written  by  J.  R.  Powell.   This  program,  written  for 
Whirlwind,  is  a  Fourier  solution  for  critical  mass  of  a  highly 
enriched,  homogeneous,  heavy  water  moderated  reactor;  methods 
were  developed  to  adapt  it  to  the  low  enrichment,  heterogeneous, 
light  water  case,  and  to  compute  the  necessary  inputs. 

3.  Computation  of  the  factor  by  which  the  critical  en- 
richment of  a  given  reactor  would  have  to  be  increased  in  order 
to  attain  a  specified  core  life. 

A  correlation  was  obtained  between  the  adjusted  Powell 
Fourier  program  and  experimental  data;  the  methods  selected 
for  calculation  of  resonance  escape  probability,  fast  fission 
factor,  and  disadvantage  factor  for  use  in  the  program  were 
also  confirmed  by  comparison  of  results  with  experimental  data. 

Ranges  of  parametric  variation  were  chosen  as  follows; 

1.  Fuel  elements  to  be  rods,  0.2,  0.4,  and  0.6  Inches 
in  diameter. 

2.  Ratio  of  volume  of  moderator  to  volume  of  fuel  to 
be  1.5,  2.0,  3.0,  4.0,  5.0,  and  6.0. 

3.  Cladding  to  be  of  stainless  steel  and  of  zirconium, 
0.015  inches  thick. 

■k.   Fuel  to  be  uranium  oxide,  UO2;  quantity  and  enrich- 
ment to  be  determined. 

5.  Cores  to  be  right  circular  cylinders  not  to  exceed 
13  feet  in  diameter  and  core  heights  to  be  4,  6,  and  8  feet. 

6.  Cores  to  be  surrounded  by  an  effectively  Infinite 
water  reflector. 

7.  Maximum  fuel  temperatures  4800°F.  cmor* 
Coolant  pressure;  2000  psia.  Inlet  temperature;   ^°7  F 

Maximum  t  emperature;  623  °F 
Average  reactor  temperature;  508°F 

8.  Thermal  power  output;  400  megawatts. 

For  each  case,  an  appropriate  core  diameter  and  total 
quantity  of  uranium  were  computed  from  the  heat  transfer  charac- 
teristics associated  with  the  parameters.   With  this,  resonance 
escape  probability,  fast  fission  factor,  and  cross  sections 
modified  for  disadvantage  factor,  the  machine  calculation  was 
entered.   The  output  of  the  machine  calculation  was  critical 
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mass  of  U-235;  this  was  compared  with  the  total  quantity  of 
uranium  previously  computed  to  give  the  required  critical  en- 
richment. 

Analysis  of  the  uranium  quantity  and  enrichment  data 
thus  obtained  indicated  that? 

1.  Total  quantity  of  uranium  required,  a  function  of 
the  heat  transfer  characteristics  of  the  unit  cell,  is  ap- 
proximately proportional  to  the  square  of  the  rod  diameter; 
the  allowable  heat  generation  rate  per  unit  length  of  the  cen- 
tral rod  is  very  nearly  the  same  regardless  of  the  rod  diam- 
eter, and  hence  the  number  of  rods  required  is  very  nearly 
the  same  for  all  rod  diameters  at  a  given  height. 

2.  Critical  enrichment  is  principally  a  function  of  the 
unit  cell  parameters  and  is  nearly  independent  of  the  core 
size  and  shape  between  diameter  to  height  of  core  ratios  of 
0.8  and  3. 

3.  For  each  rod  diameter  and  cladding  material  an  opti- 
mum moderator  to  fuel  volume  ratio  exists  at  which  the  criti- 
cal enrichment  is  a  minimum.  For  the  diameters  tested,  these 
were  as  follows? 


Rod  Diameter 


Optimum  Vi/v0 
Stainless  Steel  clad  Zirconium  clad 


3,0 


0.6  inch        3-0 
0.i+  inch        3-75 
0.2  inch        k.$ 

k.   An  optimum  rod  diameter  at  which  minimum  critical 
enrichment  is  lowest  exists  for  rods  clad  with  each  cladding 
material.   For  zirconium  clad  rods,  this  occurs  at  about  0.&5 
inches  rod  diameter  for  stainless  steel  clad,  at  slightly 
over  0.6  Inches. 

Rod  Diameter        Minimum  critical  enrichment 

Stainless  Steel  clad  Zirconium  clad 

11 1   1 1 1 1  1      11  1  1    » i«  1  ■  1  ■  1  1  ■  ■  »■■! 

0.6  inch         „99  .89 

O.k   inch        1.12  .85 

0.2  inch        1.5  1.13 

5.  Critical  enrichments  are  lower  for  zirconium  clad 
rods  than  for  stainless  steel  clad,  of  the  order  of  \%  en- 
richment compared  to  the  order  of  2%. 

6.  Optimum  critical  enrichment  does  not  necessarily  in- 
dicate the  optimum  core;  other  factors  which  must  be  considered 
are  total  quantity  of  uranium,  cost  data  for  uranium,  cladding 
materials,  and  fabrication  of  fuel  elements  and  pressure  ves- 
sels, and  the  additional  enrichment  necessary  to  attain  a 
specified  core  life. 
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NOMENCLATURE 
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A.  General 


Symbol 

A 

A 

a 

B2 

c 
D 

E 
E 

F 
F 


H 
h 

I 

J 
K 


Description 

Atomic  weight 

Coolant  channel  cross  sectional  area 

Rod  heat  transfer  area 

Buckling-a  measure  of  the  bending 
of  the  neutron  flux  at  a  point  In 
a  critical  reactor 

Coolant  specific  heat 

Diffusion  coefflclent-the  proportlon- 
allty  factor  between  neutron  current 
and  neutron  flux  gradient 

Energy 

Dlmensionless  coefficient-related  to 
the  computation  of  the  disadvantage 
factor.   A  function  of  ^ and  R.   See 
Appendix  A. 


Hot  channel  factor-heat  transfer 
safety  factor 

Thermal  utlllzatlon-the  proportion  of 
thermal  neutrons  absorbed  in  fuel  to 
the  total  absorbed 

Core  height 

Film  heat  transfer  coefficient 

Modified  Bessel  function  of  the  first 
kind 

Bessel  function  of  the  first  kind 

Modified  Bessel  function  of  the  second 
kind 


Dimensions 


ft' 


ft 


cm 


-2 


BTU/lb-F 


Mev 


cm 


BTU/hr-ft2-F 


iv 


M' 


N 


Infinite  multiplication  factor-the 
ratio  of  the  average  number  of  neutrons 
produced  in  each  generation  to  the 
average  number  of  corresponding  neu- 
trons absorbed 

Thermal  conductivity 

(Diffusion  length)  -l/6th  of  the  mean 
square  distance  traveled  by  a  thermal 
neutron  before  being  absorbed 

Migration  area°l/6th  of  the  mean 
square  distance  from  birth  of  fission 
neutron  to  capture  of  a  thermal  neutron 

6.025  x  1023yQ 


BTU/hr-ft-°F 
cm2 

2 
cm 

atoms/  cm^ 


Q 

q 

q" 

R,r 

t 

U 

V 

x 


Dlmenslonle ss  coef f icl ent°f unction 
related  to  temperature  differences  in 
heat  transfer.   See  appendix  G. 


Resonance  escape 


the 


probability  that  a  neutron  will  not  be 
absorbed  by  U-238  in  the  resonance 
energy  region 

Dlmenslonle ss  coefflclent-f unction 
related  to  heat  transfer  resistance 

Slowing  down  density-rate  at  which 
neutrons  reach  thermal  energy  per 
unit  volume 

Heat  transfer  rate.  See  Appendix  G. 

Radius 

Temperature 

Total  heat  transfer  coefficient 

Volume 

Ratio  of  volume    of   uranium  4  clad  to 
water  in  a  unit   cell 

Fast  fission  factor-the  proportion  of 
total  fissions  at  both  fast  and  thermal 
energies   in  both  U-235  and  U-238  to  the 
number  of  thermal  fissions  in  U-235 


neutrons/ 

3 
cm-7- sec 

BTU/hr-ft2 

Op 

BTU/hr-ft2- °F 
cm  3 


X 

r 

o~k 


•tr 


1 


Reciprocal  diffusion  length-  =  l/L 

Reflector  savlngs-the  decrease  in 
critical  size  of  a  reactor  due  to 
a  reflector 

Density 

Microscopic  absorption  cross  sectlon- 
the  probability  that  one  atom  of 
material  will  absorb  a  neutron 

Microscopic  fission  cross  section- 
the  probability  that  one  atom  of 
material  will  be  fissioned  by  a 
neutron 

Microscopic  scattering  cross  section- 
the  probability  that  one  atom*  of 
material  will  scatter  a  neutron 

Microscopic  transport  cross  section- 


the  product  of 


or 


and  (l-  jfQ) 


Macroscopic  absorption  cross  section- 
the  product  of  the  number  of  atoms  of 
material  per  unit  volume  and  cgt or  that 
material 

Macroscopic  fission  cross  section- 
the  product  of  the  number  of  atoms 
of  material  per  unit  volume  and  Cp 
for  that  material  » 

Macroscopic  transport  cross  section- 
the  product  of  the  number  of  atoms  of 
material  per  unit  volume  and Gj£ for 
that  material 

Age-l/6th  of  the  mean  square  distance 
traveled  by  a  neutron  in  being  slowed 
down  from  fission  to  thermal  energy 

Average  number  of  fast  neutrons  re- 
leased per  slow  neutron  fission 

Neutron  flux-the  product  of  the 
number  of  neutrons  per  unit  volume 
and  their  velocity 


cm 
cm 

gm/cm^ 
barns 

barns 


barns 

barns 
cm"1 

cm_l 


-1 


cm 


cm2 


neutrons/ 
cm2- sec 
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B.  Specific 

Symbol  Description 

»0        rod  heat  transfer  area,  midpoint  to  end  of  core 

af        rod  heat  transfer  area  midpoint  extended  through 
reflector  savings 

k^  cross  sectional  area  of  unit  cell 

Fq  heat  generation  rate  hot  channel  factor 

F£«y  temperature  rise  hot  channel  factor 

F  q  heat  transfer  coefficient  hot  channel  factor 

kc  clad  thermal  conductivity 

kf  fuel  thermal  conductivity 

P^  reactor  thermal  power 

q        average  power  developed  per  rod  averaged  over  core 

q0        average  power  developed  in  a  rod  averaged  longitudi- 
nally only 

qjj        heat  transfer  rate  at  midpoint  of  the  central  rod 

q"        heat  transfer  rate  at  midpoint  of  the  central  rod 


mBO 


required  to  produce  burnout  of  clad 


qM  heat  transfer  rate  at  midpoint  of  the  central  rod 

mLB  required  to  produce  local  boiling  at  outlet 

R0  radius  of  fuel  In  a  unit  cell  or  to  inside  of  cladding 

Rl  outside  radius  of  unit  cell 

Rc  radius  to  outside  of  cladding 

R  core  radius 

R'  core  radius  plus  reflector  savings 


vii 


A  tc      temperature  rise  of  coolant  in  core  from  inlet  to 
maximum  temperature 

Z\  tf      temperature  rise  from  coolant  inlet  to  centerllne 
of  fuel  rod 

Atwi     temperature  rise  from  coolant  inlet  to  cladding- 
fuel  interface 

/\  two     temperature  rise  from  coolant  inlet  to  outside  of 
cladding 

Uc       effective  heat  transfer  coefficient  for  bond,  clad, 
and  scale 

Vo       volume  of  fuel  in  unit  cell 

V1       volume  of  clad  4-  moderator  in  unit  cell 

V        volume  of  clad  in  unit  cell 
c 

^w,H20   volume  of  water  in  unit  cell 
Vu       volume  of  uranium  in  unit  cell 
Vs       volume  of  structure  and  clad  in  unit  cell 
~)(0  fuel  reciprocal  diffusion  length 

"yt  moderator  reciprocal  diffusion  length 

-g  fuel  macroscopic  absorption  cross  section 

<""  -.      moderator  macroscopic  absorption  cross  section 
fuel  macroscopic  transport  cross  section 
moderator  macroscopic  transport  cross  section 
slowing  down  cross  section 
^,       Fermi  age 
^*2       diffusion  age 
<P  0       volume  average  slow  neutron  flux  in  fel 

volume  average  slow  neutron  flux  in  moderator 
volume  average  slow  neutron  flux  in  clad 
disadvantage  factor 


.otr 

•  II 


ft/fl 


I. 

INTRODUCTION 

The  category  of  nuclear  reactors  known  as  pressurized 
water  reactors  has  emerged  from  the  multitude  of  possible  types 
of  power  producing  reactors  as  one  of  the  most  promising  from 
a  feasibility  viewpoint.  It  appears  that  over  the  next  sev- 
eral years  it  will  continue  to  be  an  important  type  both  for 
ship  propulsion  and  for  shore  power  station  use,  and  hence  a 
study  of  the  effects  of  varying  some  of  the  basic  core  para- 
meters was  considered  to  be  a  fruitful  undertaking. 

A  pressurized  water  reactor  is  distinguished  by  the  fol- 
lowing characteristics  ? 

1.  Water  is  used  as  both  the  coolant  and  the  moderator. 

2.  It  is  heterogeneous;  that  is,  it  is  made  up  of 
segregated  regions  of  fuel  and  coolant-moderator, 
separated  by  a  corrosion  resistant  M cladding. M 

3.  The  water  is  kept  under  high  pressure  to  prevent 
boiling. 

4.  The  great  majority  of  fissions  is  produced  by 
neutrons  at  thermal  energy ,  that  is,  neutrons 
which  are  in  thermal  equilibrium  with  the  atoms 
of  the  surrounding  medium. 

The  first  full  scale  power  reactor  to  operate,  the  pro- 
totype of  the  propulsion  reactor  in  the  U.  S.  submarine 
Nautilus,  was  of  this  type.   Since  then,  in  addition  to  the 
successful  plant  actually  installed  in  the  Nautilus,  the 
pressurized  water  concept  has  been  accepted  as  standard  for 
U.  3.  submarines  and  thirteen  other  boats  are  building,  in 


the  design  state,  or  authorized.   Aircraft  carrier  and  cruis- 
er plants  are  also  being  planned.   The  pressurized  water 
reactor  has  also  taken  a  leading  position  in  the  U.  S.  for 
shore  power  station  uses   the  first  civilian  power  reactor 
to  operate  will  be  the  PWR  at  Shipping  port,  Pennsylvania, 
and  among  the  ones  to  follow  it  will  be  the  Yankee  Atomic 
Electric  Company's  design,  to  be  built  at  Rowe,  Massachusetts. 

Similar  reactors  are  being  built  in  other  countries;  a  list 

1. 

of  those  now  planned  follows;   ' 

Army  Package  Power  Reactor 
Belgian  Thermal  Reactor 
Dominican  Republic  Reactor 
Consolidated  Edison  Reactor 
University  of  Florida  Reactor  U 

The  uranium  fuel  in  all  these  designs  is  formed  into 
elements,  such  as  plates  or  rods,  which  are  small  compared  to 
the  size  of  the  reactor  as  a  whole  but  very  large  compared  to 
atomic  distances  and  of  the  same  order  of  magnitude  as  the 
mean  free  path  of  a  neutron  in  uranium,  about  half  an  inch. 
These  elements  are  regularly  spaced  in  a  vessel  through  which 
water  can  flow,  filling  the  spaces  between  the  elements.   Each 
element  with  its  share  of  the  surrounding  water  can  be  thought 
of  as  a  cell,  and  the  whole  reactor  as  an  aggregate  of  such 
cells,  or  a  cell  lattice.   Some  of  the  nuclear  and  thermal 
properties  of  the  reactor  are  properties  of  the  cell  parameters 


Fuel 

Heat 
Power 

Critical  Date 

TJ02 

10  mw 

Early  1957 

U 

260 

1957 

vo2 
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u 

500 

July  1959 

u 

10 

December  1959 

1." Nuclear  Reactor  Data  Handbook  No.  2",  The  Raytheon  Mfg.  Co. 
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while  some  are  properties  of  the  lattice.   The  properties  of 
basic  interest  ares—1- 

1.  Nuclear  properties  of  the  unit  cell. 

2 

a.  The  infinite  multiplication  factor—1-  k00.   In 

order  to  maintain  steady  state  operation,  the 
number  of  neutrons  in  successive  generations 
must  be  conserved.   That  is? 

Production  =  Absorption  +  Leakage 
k00  is  the  number  of  thermal  neutrons  produced 
per  thermal  neutron  absorbed  if  no  leakage 
takes  place.   It  can  be  thought  of  as  the 
factor  relating  the  numbers  of  neutrons  in 
successive  generations  in  a  medium  consisting 
of  an  infinite  repetition  of  the  unit  cell. 
It  is  generally  recognized  as  the  product  of 
four  factors,  which,  with  the  energy  cycle 
followed  by  the  neutrons,  are  explained  below. 
V     '   The  number  of  fast  neutrons  produced 
per  thermal  neutron  absorbed  in  the 
fuel.   High  speed  neutrons  with  a  mean 
energy  of  about  1  mev  are  produced  by 
absorption  of  thermal  neutrons,  those 
which  have  been  so  slowed  down  by 


T~.      fhese  properties  are  described  in  more  detail  in  various 
references,  and,  in  some  cases,  in  the  Appendices  to  this 
thesis.   These  references  are  cited  in  each  case. 

2.   See  S.  Glasstone  and  M.C.  Edlund,  MThe  Elements  of  Nuolear 
Reactor  Theory*,  Van  Nostrand  (1952) (hereafter  cited  as 
"Glasstone  and  Edlund"),  Section  h.57   ff» 


4. 

successive  colllsons  that  their  energies 
depend  only  on  the  temperature  of  the 
surrounding  medium.   Not  all  those  absorbed 
in  the  fuel  produce  fissions;  those  that 
do,  produce  on  the  average  2.^6  fast  neu- 
trons and  about  192  mev  of  energy. 
€        :   The  fast  fission  factor.   Neutrons  with 
energies  greater  than  about  1  mev  are 
capable  of  producing  a  small  number  of 
fissions  in  U-238.  £  is  the  ratio  of  the 
total  rate  of  neutron  production  to  the 
production  rate  from  thermal  fission  of 

U-235. 
p    :    The  resonance  escape  probability.   U-238 

has  a  high  probability  of  absorption  with- 
out fission  for  neutrons  in  the  "resonance11 
energy  range,  from  about  10  to  100  ev.   p 

is  the  probability  that  a  fission  neutron 

in  slowing  down,  will  not  be  so  captured. 

f   :   The  thermal  utilization  is  the  proportion 

of  all  the  thermal  neutrons  absorbed 
which  are  absorbed  In  the  fuel. 

b.  The  disadvantage  factor,—1"  /tf)0  •   Thermal 

neutrons  are  produced  in  the  moderator  region 
as  fast  neutrons  resulting  from  fissions  in 


1.   See  Glasstone  and  Edlund,  9*^8,  and  Appendix  A.   Also  see 
A.  L.  Kaplan,  "Theoretical  Studies  of  Neutron  Flux  Distri- 
bution Expected  in  the  MIT  Nuclear  Research  Reactor",  MIT 
S.M.  thesis,  1955. 


the  fuel  region  are  slowed  down  "by  elastic 
collisions  with  the  hydrogen  nuclei. 

The  neutron  flux,—1-  d   ,  expresses  the  length 
of  neutron  paths  covered  per  unit  time  per  unit 
volume.   This  depends  both  on  the  number  of 
neutrons  and  on  their  velocity.   This  is  shown 
in  the  units  of  flux: 

neutrons  x  -»  x  X         1 


sec*   cm3     cm2~sec 
The  number  of  interactions  taking  place  per 
unit  time  and  volume  is  the  product  of  the 
length  of  neutron  paths  and  the  probability 
of  interaction  per  unit  path  length,^. 

no.  of  Interactions    -  ^    $ 
unit  time-unit  volume         » 

^  in  turn  consists  of  the  number  of  nuclei 
per  unit  volume ,  N,  times  the  reaction  proba- 
bility per  nucleus pC?-".     Thus, 

no.  of  interactions    *  ^(fi    *  fj&-*  S 
unit  time-unit  volume 

The  fast  flux  produced  by  fissions  in  the 

fuel  is  slowed  down  in  the  moderator,  producing 

thermal  flux.   As  the  thermal  flux  thus  produced 

reenters  the  fuel,  it  is  continuously  reduced 

as  the  center  of  the  fuel  element  is  approached 

since  ^  for  absorption  is  greater  in  the  fuel 

than  in  the  moderator  and  since  essentially 

no  slowing  down  takes  place  in  the  fuel  to 


1.  See  (Jlasstone  and  Edlund,  3.49. 


replenish  the  supply  of  thermal  neutrons.   The 
general  pattern  of  thermal  flux  In  a  oell  Is 
shown  in  Figure  1. 

The  fuel  near  the  center  of  the  element , 
thus s  meets  a  lower  flux  than  the  average  over 
the  whole  reactor,  with  a  consequent  reduction 
of  reaction  probability  below  that  which  would 
be  exhibited  by  a  homogeneous  mixture  of  the 
same  materials  in  the  same  average  thermal 
flux.   The  ratio  of  the  average  thermal  neutron 
flux  in  the  moderator  to  that  in  the  fuel  is 
the  disadvantage  factor. 


c.  Fermi  Age,   3f  .   Since  the  slowing  down  process 
is  a  complex  one,  involving  fission  neutrons 
with  a  distribution  of  Initial  energies  taking 
part  in  a  large  number  of  elastic  collisions 
in  which  their  energy  loss  varies  (for  a  hydro- 
gen containing  medium)  between  zero  and  100#, 
it  is  represented  by  an  experimentally  deter- 
mined number  which  is  characteristic  of  the 


1.  See  Glasstone  and  Edlund,  6.117  to  6.145,  and  Appendix  F. 


proportions  of  the  various  materials  in  the 
reactor;  since  these  materials  and  proportions 
are  repeated  in  all  the  cells,  it  can  thus  "be 
considered  a  property  of  the  cell.   Physically, 
the  Fermi  age  is  one  sixth  of  the  average 
squared  crow=f light  distance  travelled  by  a 
neutron  during  the  process  of  slowing  down 
from  fission  to  thermal  energy. 
2.  Nuclear  properties  of  the  reactor  as  a  whole. 

a.  Flux  distribution,—1"  &[+).      Due  to  the  effects 


of  leakage  and  changes  of  composition  across 
the  lattice,  the  flux  varies  over  each  dimen- 
sion of  the  reactor,  being  high  near  the  center 
and  in  regions  of  high  multiplication  factor, 
and  low  near  the  edges  and  in  regions  of  low 
multiplication  factor.   It  drops  to  zero  a 

small  distance  outside  the  physical  boundaries 

2. 

of  the  reactor  at  the  "extrapolated  boundary. "— 1- 

The  functions  which  express  these  patterns  are 

as  follows? 

Reactor  Shape      Dimension       Function 

Cube         parallel  to  an  edge   003/77"  x  \ 
of  length  H  ^2  Wz) 


Right  clrcu-  radial  (radius  s  R)  J0  (2.-  Wr J 

lar  cylinder  \  R   / 

axial  (height  =  H)  cos/ff  x_\ 

Sphere        radial  (radius  sR)  TTr  slnhrv) 


1.  See  Qlasstone  and  Edlund,  7.36  to  7«6l 

2.  See  Glass tone  and  Edlund,  5.40 
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A  reflector  is  frequently  provided  to  reflect 
neutrons  back  into  the  core  which  would  other- 
wise leak  from  it.   When  this  is  done,  the 
average  flux  in  the  active  core  is  increased, 
as  shown  in  Figure  2,  and  the  flux  pattern 
is  extended  to  the  reflector  boundary,  c. 
The  flux  in  the  reflector  has  a  different 
curvature  from  that  in  the  core;  but  the 
core  flux  can  be  approximated  by  extending 
the  point  at  which  its  characteristic  function 
becomes  zero  by  a  distance  such  as  ab,  called 
the  reflector  savings.—^ 

fiBFLECTED 
JZLLECTj 


Figure  2 

2. 

b.  Buckling,— ~  b2o  is  a  measure  of  the  amount  of 

curvature  of  the  neutron  flux  pattern.   The 
buckling  which  is  required  to  make  a  given 
reactor  critical  can  be  found  from  the  shape 
and  dimensions  of  the  reactor ,  and  that  which 
can  be  obtained  by  a  given  combination  of 
materials  can  be  found  from  the  properties  of 
the  materials.   The  condition  for  criticallty 


1.  See  Glasstone  and  Edlund,  8.25,  and  Appendix  F. 

2.  See  Glasstone  and  Edlund,  7»20  to  ?.60. 
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of  the  reactor  Is  that  these  two  "Ducklings  be 
equal.   Buckling  can  also  be  used  to  find  the 
proportion  of  neutron  leakage  from  a  given 
reactor .   The  probability  of  non-leakage  dur- 
ing the  slowing  down  process  is  e=B^(  and  the 
ratio  of  thermal  leakage  to  thermal  absorption 

is  L2B2,  where  L  is  the  thermal  diffusion 

1. 
length. — " 


k 
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Figure  3 
3»   Thermal  properties  of  the  unit  cell. 

a.  Maximum  heat  generation  rate-—-  q" 


The  rate 


at  which  heat  can  be  generated  in  a  fuel  ele- 
ment is  limited  by  Its  heat  transfer  proper- 
ties.  Along  an  element,  the  heat  generation 
rate  is  proportional  to  the  fission  rate  which 
is  proportional  to  the  thermal  neutron  flux, 
which  will  vary  from  zero  at  the  ends  of  the 
reflector  savings  to  a  maximum  halfway  between 
them.   For  equal  reflector  savings  at  each 


1.  See  Glasstone  and  Edlund,  5*62  to  5»91 

2.  See  W.  M.  Rohsenow,  J.  Lewlns,  and  J.  P.  Barger,  "Steady 
State  Temperature  Distribution  in  a  Nuclear  Reactor  with 
End  and  Center  Fed  Coolant",  NT-269,  MIT,  hereafter  cited 

as  Rohsenow,  Lewlns,  and  Barger,  pages  1=4.  Also  see  Appendix  Q-. 
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end  of  the  elements ,  this  maximum  will  also 
be  at  the  geometrical  center  of  the  active 
core.   The  coolant  flows  along  the  element 
and  its  temperature,  t0,  rises  as  it  does 
so  to  a  maximum  at  the  outlet,  while  the 
temperature  at  the  element  centerline,  t£  , 
rises  to  a  maximum  at  a  point  between  the 
middle  of  the  element's  length  and  the  out- 
let, then  drops  off  due  to  the  reduction  in 
the  heat  generation  rate.   tie  rate  at  which 
heat  can  be  generated  in  a  given  element  is 
limited  by? 

(1)  the  maximum  fuel  centerline  temperature 
must  be  maintained  below  the  melting 
point  of  the  fuel  material,  and,  for 
uranium  metal  fuel,  below  the  phase 
change  at  660  °F. 

(2)  the  maximum  coolant  temperature  must 
be  less  than  the  boiling  point  of  the 
coolant.   Boiling  will  change  the  water 
density,  changing  the  density  of  hydro- 
gen atoms  in  the  moderator  and  thereby 
the  slowing  down  characteristics  and 
the  reactivity  of  the  reactor.   In  ad- 
dition, boiling  will  reduce  the  heat 
transfer  rate  away  from  the  element, 
leading  to  a  temperature  rise  which  may 
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be  sufficient  to  melt  it. 
(3)  The  temperature  of  the  cladding  will 
follow  a  pattern  between  those  of  the 
element  centerllne  and  the  coolant. 
This  must  also  be  maintained  below  the 
melting  point  of  the  cladding. 
To  improve  the  thermal  efficiency  of  the  plant, 
it  is  desirable  to  make  the  coolant  outlet 
temperature  as  high  as  possible.   Since  this 
is  limited  by  the  boiling  point  of  the  coolant, 
the  allowable  temperature  may  be  increased  by 
increasing  the  water  pressure.   As  can  be  seen 
from  the  curve  of  Figure  4~" "  pressures  greater 
than  2000  psia  give  very  little  gain  in  tempera- 
ture; moreover;,  the  critical  point  of  water, 
at  which  the  distinction  between  the  liquid 
and  vapor  phases  is  lost,  is  at  3206. 2°F,  which 
places  an  absolute  upper  limit  on  the  pressure. 
A  more  practical  limit  is  imposed  by  the  dif- 
ficulties of  manufacturing  pressure  vessels 
of  the  required  sizes  for  pressures  greater 
than  around  2000  psia. 

7o& 
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Figure  k 


T7j7  H.    Keenan  and  F.   (J.    Keyes,    "Thermodynamic  Properties  of 
Steam",    Wiley,    (1936),    page   39- 
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4.  Thermal  properties  of  the  whole  reactor. 

a.  Overall  heat  generation  rate.—    The  heat  genera- 
tion rate  in  each  element  will  be  proportional 

to  its  q^s  which;,  in  turn,  is  proportional  to 
the  flux  distribution  in  a  plane  perpendicular 
to  the  axis  o,f  the  element .   Thus  the  maximum 
heat  generation  rate  allowed  for  an  element 
can  exist  only  in  one  element,  at  the  point 
of  maximum  flux,  and  all  other  elements  will 
operate  at  heat  generation  rates  less  than  the 
maximum.   The  ratio  between  the  maximum  and 
average  heat  generation  rates  and  the  maximum 
and  average  thermal  neutron  fluxes  will  be  the 
same . 

b.  Temperatures ,  t.   The  coolant  inlet  temperature 
is  determined  by  the  power  system  outside  the 
reactor;  the  coolant  outlet  temperature  is  the 
temperature  of  the  mixed  outlets  of  all  the 
channels »   The  average  reactor  and  coolant 
temperatures  can  be  found  from  the  variations 
in  point  temperatures  through  the  core. 

A  uranium  fuelled,  light  water  cooled  and  moderated,  pres- 
surized water  reactor  can  be  described  by  several  basic  para- 
meters: 

1.  Describing  the  unit  cells 
a.  Element  shape  and  size. 


1.  See  Appendix  G. 
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b.  Volume  ratio  of  moderator  to  fuel. 

c.  Material  and  thickness  of  cladding. 

d.  Fuel  material  and  degree  of  enrichment  In 

U-235. 
2.  Describing  the  reactor  as  a  whole: 

a.  Core  shape  and  dimensions. 

b.  Reflector  dimensions. 

c.  Temperature  limitations. 

d.  Power  output. 

Changes  in  these  parameters  act  and  interact  in  complex 
ways  to  affect  the  performance  of  the  reactor.   The  number  of 
parameters  involved  and  the  ranges  of  possible  variation  make 
experimental  determination  of  the  optimum  combination  for  given 
requirements  a  long,  tedious,  and  expensive  process. 

Exponential  experiment s—1  have  been  performed  by  H.  J. 
Kouts  at  the  Brookhaven  National  Laboratory  on  a  number  of 
the  combinations—  and  some  of  these  have  been  continued  to 
criticality  by  S.  Krasik  and  A.  Radkowsky  (of  the  U.  S.  Atomic 
Energy  Commission)  at  the  Westinghouse  Electric  Company's  Bettis 
Plant. ^-    Both  these  series  involved  uranium  metal  fuel  elements; 


1.  See  Glasstone  and  Edlund,  9. 90  to  9. 111. 

2.  Herbert  Kouts,  "Intracell  Flux  Traverses  and  Thermal  Utiliza- 
tions for  1.15$  Enriched  Uranium  Rods  in  Ordinary  Water* , 
Brookhaven  National  Laboratory  Report  No.  1987,  August  11, 
195^.   H.  J.  Kouts,  J.  Chernick,  and  I.  Kaplan,  "Exponential 
Experiments  on  Light  Water  Moderated  1  per  cent  U-235  Lat- 
tices", Brookhaven  National  Laboratory  Report  No.  209^, 
November  28,  1952. 

3.  S.  Krasik  and  A.  Radkowsky,  "Pressurized  Water  Reactor 
Critical  Experiments",  paper  no.  601,  presented  at  the 
Geneva  Conference  on  the  Peaceful  Uses  of  Atomic  Energy, 
1955. 
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a  similar,  though  shorter,  series  of  critical  experiments  on 
uranium  oxide  elements  has  been  made  by  J.  R.  Brown,  A.  Z. 
Krantz,  and  others  at  Bettis~ 

These  data  are  not  by  themselves  sufficient  to  select 
an  optimum  configuration  for  a  given  set  of  requirements,  but 
they  are  sufficient  to  provide  experimental  verification  of 
calculations. 

The  experimental  results  include  critical  buddings, 
and,  in  those  continued  to  criticality,  critical  number  of 
rods,  disadvantage  factors,  fast  fission  factors,  resonance 
escape  probabilities,  and  thermal  utilizations  for  a  series 
of  rod  diameters  and  volume  ratios  at  a  limited  number  of  en- 
richments between  1  and  1.3#»   In  addition,  Routs"  work  in= 
eludes  some  extremely  valuable  plots  of  relative  neutron  flux 
along  various  lines  through  the  cells.   All  the  fuel  was  in 
rod  form  and  clad  with  aluminum,  with  active  core  heights  of 
48  Inches.   All  cores  were  surrounded  by  effectively  infinite 
water  reflectors. 

For  the  purpose  of  theoretical  investigation  of  a  number 
of  different  configurations,  several  mathematical  models  are 
available,  varying  in  both  complexity  and  accuracy.   The 


1.  J.  R.  Brown,  B.  H.  Noordoff,  F.  S.  Frantz,  J.  J.  Volpe,  and 
D.  R.  Harris,  "Measurement  on  Low  Enrichments  Light  Water 
Moderated  Critical  Assemblies  of  U  and  0*02%  paper  no. 
XVII-5  presented  at  the  Chicago  Meeting  of  the  American 
Nuclear  Society,  1956. 

A.  Z.  Krantz,  G.  0.  Smith,  D.  Klein,  and  W.  Baer,  Mf,  p, 
and   Measurements  on  Slightly  Enriched  Light  Water  Moder- 
ated U02  Lattices",  paper  no.  XVII-6  presented  at  the 
Chicago  Meeting  of  the  American  Nuclear  Society,  1956. 
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simplest  of  these,  the  "one  group"  equatlon-7'  assumes  that  all 
neutrons  are  at  thermal  energy.   The  critical  condition  is 


1  = 


1  +  B2L2 


This  is  relatively  simple  mathematically,  but  poor  in  accuracy. 

Somewhat  better  in  accuracy  and  not  much  more  difficult  in 

2 
computation  is  the  "modified  one  group"  theory^-1-  which  in- 
cludes a  factor  for  leakage  during  the  slowing  down  process. 
The  critical  condition  may  be  expressed  as 

K>     e  -#* 

1  4  bV 
the  Fermi  Age  variant,  or 


1  = 


1   t  B2M2 

where  M2,  the  migration  area, 2-  is  (L  ^(T)i    for  a  large 
reactor. 

More  complex  is  the  "two  group"  theoryj""  which  assumes 
a  fast  and  a  thermal  group  of  neutrons  with  average  reactor 
properties  for  each.   Four  differential  equations  are  formed, 
expressing  fast  and  thermal  flux  in  the  core  and  the  reflector 


1.  See  Glasstone  and  Edlund,  8.15  to  8.37. 

2.  See  Glasstone  and  Edlund,  7.20  to  7.27,  and  7.62  to  7.64 

3.  See  Glasstone  and  Edlund,  7^63 

4.  See  Glasstone  and  Edlund,  8.38  to  8.58 
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Where  the   subscripts  1  and  2  refer  to  the  fast  and   thermal 
fluxes  respectively  and  the   subscripts  c  and  r  refer  to  the 
core  and  reflector  respectively. 

These  are   solved  by  means   of  boundary  conditions  which 
require  both  fluxes  to  be   continuous   across  the  core-reflector 
boundaries  and  both  fluxes  to  go   to  zero  at  the   outer  reflec- 
tor extrapolated  boundary.      These   solutions  lead  to  a  critical 
determinant. 
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D  =  diffusion  coefficient1 

V  s  Fermi  Age 

Coefficients  for  other  geometries  are  given  in  G-lasstone  and 

Edlund,  sections  8.^6  and  8.51* 

This  treatment  is  capable  of  considering  only  one  reflec- 
tor. But  in  practical  reactors,  reflectors  are  usually  pro- 
vided for  all  exposed  faces.   This  situation  can  be  handled 
with  moderate  accuracy  by  Hill°s  Approximation—  in  which  the 
reactor  is  considered  separately  for  each  reflected  dimension 
and  the  solutions  superimposed.   In  this  approximation,  the 
corners  are  not  properly  considered,  but  the  flux  is  generally 
low  in  those  areas  and  the  resulting  error  not  excessive. 

It  is  possible  to  reduce  the  labor  Involved  in  critical 
computations  considerably  by  use  of  a  digital  computer.   These 
calculations  are  particularly  useful  in  an  investigation  such 
as  this,  where  a  large  number  of  similar  problems  will  have 
to  be  solved  by  identical  methods,  for  once  the  program  is 
written,  it  cart  be  used  to  solve  any  number  of  cases. 

Two  such  programs  for  Whirlwind,  the  high  speed  digital 

computer  at  MIT  were  available,  both  written  by  James  R. 

2 
Powell  as  part  of  his  MIT  thesis  for  the  Sc.D.  degree.— - 

1.  G-lasstone  and  Edlund  5«7  and  5.20  to  5.2^ 

2.  Powell,  James  R. ,  "Nuclear  Characteristics  of  Heavy  Water 
Moderated,  Enriched,  Homogeneous  Reactors",  Sc.D.  thesis, 
MIT,  1957 
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These  programs  were  based  on  the  assumptions  of  a  heavy  water 
moderator,  highly  enriched  uranium,  and  a  homogeneous  core. 
One  of  these  programs  was  a  computer  representation  of  the  two 
group  Hill's  approximation  method  already  discussed;  the  other 
was  based  on  a  Fourier  treatment  of  the  diffusion  equations. 

The  theory  of  the  Fourier  treatment  will  be  briefly  out- 
linedhere.— -  A  detailed  description  of  the  theory  for  the 

Fourier  program  is  contained  in  Powell's  thesis.   The  thermal 
diffusion  equation  states  that  the  losses  due  to  leakage  and 
absorption  must  equal  the  number  of  neutrons  s lowing  down  to 
thermal  energy  for  a  critical  reactor. 


2il  -  &A 


D  D  D 

Leakage   absorption   production 


(1) 


Each  of  the  terms  in  equation  (1)  can  be  expanded  in  a  Four- 
ier  series  with  elgenf unctions  obtained  from  the  wave  equa- 
tion 

V%  T   I  i 

to  satisfy  boundary  conditions.   This  thermal  flux  expansion 

/ 
would  take  the  form  of 


+=0 


ai*i 


T~.      Based  on  J.  R.  Powell's  Thesis,  Section  III 
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Since  absorption  and  slowing  down  vary  across  the  reactor  and 
depend  on  factors  other  than  flux,  two  further  expansions  are 
required 


dTFT 


+*o 


Since  ^r a(r)  i8  not  constant,  but  varies  as  a  step  function 
of  (r)  from  region  to  region  an  additional  equation  is  writ- 
ten relating  the  ^  with  the  regions.   The  coefficients  of 
these  expansions  are  then  related  to  the  Fourier  coeffici- 
ents of  the  thermal  flux  expansion.   Using  these  relations, 
an  infinite  series  of  equations  is  obtained  with  one  equation 
for  each  unknown  Fourier  coefficient.   For  non-trivial  solu- 
tions the  determinant  of  the  coefficients  must  equal  zero. 
This  constitutes  the  critical  equation.   In  the  program  the 
critical  mass  is  found  using  a  limiting  process,  by  finding 
the  mass  of  U-235  needed  to  make  the  determinant  of  the  first 
two  rows  and  columns  vanish,  then  the  determinant  of  the  first 
three  rows  and  columns  vanish  and  noting  that  the  mass  approaches 
a  limiting  value  as  the  number  of  rows  and  columns  Increases. 

Powell  points  out  two  major  advantages  to  the  Fourier 
solution  of  diffusion  equations.  The  first  of  these  is  the 
ability  to  represent  the  slowing  down  behavior  of  fast  neu- 
trons as  a  convolution  of  an  age  and  diffusion  kernel.  The 
Fermi  age  kernel  is  expressed  by 
Pi  =  exp(  -?"iBf) 
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and  the  diffusion  kernel  is  expressed  by 

p2=    — i _ 

i  +  r2  b±z 

By  postulating  that  the  neutrons  followed  Fermi  age  theory 
until  they  reached  energy  E1   and  then  diffused  away  from  the 
source  until  a  collision  reduced  their  energy  to  thermal 
energy,  Eth,  it  was  felt  by  J.  Powell  that  effectively  two 
fast  groups  of  neutrons  were  being  considered.   The  slowing 
down  kernel  was  then  represented  as 

P  =  t^   =  ?!  [ttl^Eil.Bi2]  P2  (EltBi2) 

2 


1  +r2  Bt' 

The  relative  size  of  ^   and  tf2   was  determined  by  matching 
the  observed  distribution  of  slowing  down  density  with  the 
calculated  distribution  so  that  no  definite  energy  boundary 
between  fast  groups  is  established. 

The  Second  major  advantage  of  the  Fourier  treatment 
pointed  out  by  Powell  was  the  relative  ease  with  which  addi- 
tional reflectors  could  be  handled.   As  a  result  of  this, 
the  Powell  program  included  the  provision  for  two  top  and 
two  bottom  reflectors  and  for  six  radial  regions,  any  one 
of  which  could  be  the  core  region.   All  cores  in  this  thesis 
are  analyzed  or  designed  on  the  basis  of  two  regions  only,  a 
core  and  an  outer  infinite  water  reflector  on  top,  bottom  and 
sides. 
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A  choice  of  two  methods  of  machine  computation  was  thus 
offered:   the  Fourier  treatment  and  the  two  group  Hill's  ap- 
proximation.  For  several  reasons,  the  Fourier  program  seemed 
preferable  for  the  purposes  of  the  thesis,  and  it  was  ulti- 
mately selected.   These  reasons  were: 

1.  The  simultaneous  solution  of  two  cases  required 
for  Hill's  approximation  makes  the  solution 
rather  lengthy.   Computer  time  required  for  this 
program  is  about  four  minutes  compared  to  two 
minutes  for  the  Fourier. 

2.  Comparisons  made  by  Powell  in  his  thesis  of  the 

calculated  results  by  both  methods  against  ex- 

perimental  data  for  highly  enriched,  heavy 

water  moderated  cores  showed  smaller  errors  In 

the  critical  masses  calculated  by  the  Fourier 

program.   Errors  in  cylindrical  geometry  varied 
from  t  .08  to  -  1k%   with  the  Fourier  and  from 

-  4.0  to  -  20#  with  the  two  group  Hill's  approxi- 
mation. 

3.  Use  of  a  two  group  program  requires  Inputs  of 
fast  group  constants.   These  are  not  well  known, 
and,  since  a  variety  of  materials  was  to  be 
used  in  the  core,  it  appeared  that  the  introduc- 
tion of  this  added  uncertainty  would  not  be  ad- 
vantageous.  The  Fourier  program,  on  the  other 
hand,  uses  only  thermal  constants  and  a  Fermi 
Age,  all  of  which  are  known  to  fair  accuracy. 
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Having  the  required  program  Inputs,  the  data  available 
for  experimental  checks,  the  practical  problems  of  fabrica- 
tion which  would  limit  the  choices  in  a  real  systen  in  mind, 
the  following  ranges  of  parameter  variation  were  selected: 

1.  In  the  unit  cell. 

a.  Fuel  elements  to  be  rods,  from  0.1  to  0.6 
inches  in  diameter. 

b.  Volume  ratios,  moderator  to  fuel,  from  1.5  to 
6.0. 

c.  Cladding  to  be  .015  inches  thick,  both  stain- 
less steel  and  zirconium  to  be  considered. 

d.  Fuel  to  be  uranium  oxide,  UO2J  quantity  and 
enrichment  to  be  determined. 

2.  For  the  whole  reactor. 

a.  Cores  to  be  right  circular  cylinders  not  to 
exceed  13  feet  in  diameter,  and  to  vary  between 
k   and  8  feet  in  height. 

b.  Cores  to  be  surrounded  by  an  18  centimeter  thick 
water  reflector  both  radial  and  top  and  bottom; 
this  is  effectively  an  Infinite  reflector. 

c.  Maximum  temperature  in  fuels  ^800°F 
Coolant  pressure:  2000  psia 
Coolant  temperatures:  Maximum:  628°F 

Inlet   :  487°F 
Average:  508°F 

d.  Thermal  power  output:  4-00  megawatts. 
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In  order  to  carry  out  the  parametric  study  utilizing 
the  machine  calculation,  a  series  of  steps  were  undertaken? 

1.  A  method  was  developed  to  convert  the  program  from 
the  high  enrichment,  homogeneous  case  with  heavy 
water  moderator  and  aluminum  structure  for  which 
it  was  written ,  to  the  case  which  was  to  "be  in- 
vestigated?  low  enrichment,  heterogeneous, 

light  water  moderated,  and  stainless  steel  or 
zirconium  structure. 

2.  A  method  was  developed  for  designing  a  consistent 
series  of  cores  for  testing  which  reflects  the 
influence  of  parameters  to  be  varied  both  on  the 
heat  transfer  and  the  nuclear  characteristics 

of  the  core,  and  supplies  the  necessary  data  for 
entering  the  machine  calculation. 

3.  A  survey  was  made  to  collect  best  values  of  the 
necessary  constants  and  basio  data,  and  to  select 
methods  for  the  various  preparatory  calculations 
required. 

4.  Correlations  were  obtained,  both  for  the  prepara- 
tory calculations  and  for  the  computer  results, 
with  the  experimental  data  available. 
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II. 

PROCEDURE 

The  use  of  the  Whirlwind  program  selected,  which  was 
written  for  a  highly  enriched,  heavy  water  moderated,  homogen- 
eous reactor,  for  computations  on  a  low  enrichment,  light  water 
moderated,  heterogeneous  reactor  involved  a  considerable  amount 
of  calculation  prior  to  entering  the  machine  computation. 

The  program  assumes  the  existence  of  only  three  materials 
in  the  core: 

1.  Fuel,  assumed  to  he  U-235. 

2.  Moderator,  assumed  to  be  D20. 

3.  An  absorber,  assumed  to  be  aluminum. 
Appropriate  values  for  the  densities  and  atomic  weights  of 
these  materials  are  built  into  the  program,.   Further,  the 
program  assumes  that  the  core  region  consists  of  a  uniform 
mixture  of  moderator  and  absorber  with  a  small  amount  of  fuel 
uniformly  distributed  in  the  mixture.   Among  the  computer 
inputs-*- '  are  values  of  the  microscopic  absorption  and  transport 
cross  sections  of  moderator  and  absorber,  the  volume  ratio  of 
absorber  to  moderator,  and  the  microscopic  absorption  and 
fission  cross  sections  of  the  fuel.   The  principal  problem 

in  utilization  of  the  program  was  to  compensate  within  these 
inputs  for  the  effects  of  heterogeneity,  the  use  of  materials 
other  than  heavy  water  and  aluminum  as  moderator  and  absorber, 
and  the  additional  presence  of  relatively  large  quantities  of 
U-238,  for  which  the  program  makes  no  provision. 

T~.    See  Appendix  E. 
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In  principle,  two  reactors,  one  heterogeneous,  the  other 
homogeneous  are  equivalent  if  the  average  values  of  macroscopic 
absorption,  fission,  and  transport  cross  sections  for  the  en- 
tire reactors  are  equal.   This  is  the  basis  of  the  derivations 
for  the  necessary  adjustments  to  the  cross  sections  as  outlined 
in  Appendix  B.   Each  value  of  a  thermal  cross  section  is  ad- 
Justed  to  have  the  same  value  as  the  effective  volume  average  value 
of  the  corresponding  cross  section  in  the  heterogeneous  reactor 
weighted  by  the  relative  flux  through  the  unit  cell.   For  ab- 
sorption in  the  fuel,  as  an  example,  this  is  expressed  as 
follows: 

^a(25)  -  N(25)^T(25)  Vo  fig      .  N'(25)c7-a«(25)  (1) 

vo  $>  +  Vc  ft 

£  a(25)  =  N(25)^(25)   vV<>    =  H(25)  YQ  <^a»  ftp  (2) 

^  a  Vo  +  V  e  _JT  V0+V  o 

% 

where  N(25)  and  o-^a(25)   are  the  number  of  atoms  per  cubic  centi- 
meter of  fuel  and  the  microscopic  cross  section  in  the  hetero- 
geneous reactor  and  N°(25)  andc7-a8(25)  are  the  number  of  atoms 
per  cubic  centimeter  of  core  and  the  microscopic  cross  section 
in  the  equivalent  homogeneous  reactor,  <7-fe°(25)  is  the  appro- 
priate computer  input. 

/fL    is  the  disadvantage  factor,  discussed  in  the 
Introduction.   Methods  of  computing  the  disadvantage  factor 
are  discussed  in  Appendix  A.   The  choice  of  method  for  use  in 
this  thesis  was  based  on  the  results  of  A.  L.  Kaplan^s  M.I.T. 


1.  See  Appendix  B. 
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thesis,  '  in  which  he  compared  the  accuracy  of  the  disadvan- 
tage factor  for  the  M.I.T.  Reactor  as  computed  from  three  in- 
creasingly sophisticated  models.   His  conclusion  was  that  the 
most  complex  of  the  three  models  offered  only  a  three  per  cent 
increase  in  accuracy  over  the  simplest,  and  that  the  added 
accuracy  did  not  justify  the  added  labor  of  computation.   The 
model  chosen,  then,  was  the  simplest  of  Kaplan °s  three;  it 
assumes  constant  fast  flux  across  the  cell  and  no  net  thermal 
leakage  from  the  cell.  * 

Since  the  program  is  based  on  pure  U-235  as  the  fuel, 
it  appeared  that  the  best  way  to  include  the  effects  of  U-238 
was  to  consider  it  as  part  of  the  absorber,  which  then  consisted 
of  both  U-238  and  the  cladding  material.   The  process  developed 
for  adjusting  the  constants  for  this  arrangement  was  as  follows s 

1.  A  macroscopic  cross  section  was  formed  for  the 
U-238  from  N(28)  in  the  fuel  andc^a' (28),  obtained 
by  a  disadvantage  factor  adjustment  as  in  equation 
(3)  above. 

2.  A  macroscopic  cross  section  was  obtained  for  the 
cladding  material.   Since  the  cladding  is  in  a 
region  of  the  unit  cell  in  which  the  flux  is  very 
nearly  the  average  value,  no  disadvantage  factor 
correction  was  included,  and  $he  macroscopic 
cross  section  consisted  simply  of  the  product 
N(clad)  x^-&(clad). 


1.  A.  L.  Kaplan,  "Theoretical  Studies  of  Neutron  Flux  Distri- 
bution Expected  in  the  MIT  Nuclear  Research  Reactor",  MIT 
S.M.  Thesis,  1955. 

2.  Also  discussed  in  Glasstone  and  Edlund,  9.47  to  9.66,  and 
in  Appendix  A. 
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3.  These  cross  sections  were  weighted  by  their  respec- 
tive volumes  and  averaged  over  the  total  volume 

of  fuel  and  cladding. 

4.  Since  the  machine  would  compute  a  macroscopic  ab- 
sorber cross  section  on  the  basis  of  the  density 
and  atomic  weight  of  aluminum,  which  were  integral 
parts  of  the  program,  the  average  macroscopic 
cross  section  obtained  in  step  (3)  was  divided 

by  N(A1). 
Mathematically,  the  process  is  described  by  the  following  equa- 
tion: 

M absorber)  -  N(28)^(28)  VQ  ±   N  ( clad  )c^(  clad)  V(clad) 
crtaosoroer;  v©  +  V(clad)  x 

AA1 


z>A1  x  6.025  x  1023 
A  step  for  spreading  the  absorber  uniformly  through  the  core 
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volume  is  part  of  the  computer  program;  a  program  input,  MxM, 
is  the  volume  ratio  of  absorber  to  moderator.   For  computing 
Mxw,  the  volume  of  fuel  rod  plus  cladding  was  used. 

A  further  adjustment  must  be  made  for  the  presence  of 
U-238  in  the  fuel  rods.   With  highly  enriched  uranium,  there 
is  neither  a  fast  fission  effect  nor  resonance  absorption^  and 
the  infinite  multiplication  factor,  k00  is  simply 

•^for  the  fuel  is  a  program  input;  f  is  computed  by  the  machine 
However,  for  low  enrichments, 

k. 


oo 


?y»f 
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In  order  to  use  the  program  for  low  enrichments,  then,  £  and 
p  were  separately  calculated,  as  described  in  Appendix  D,  and 
the  product  was  used  as  the  program  input  *fn. 

An  additional  arbitrary  factor  of  1.1  was  used  in  the 
program  input  My  M  to  bring  computed  results  into  agreement 
with  experimental  data.   The  derivation  of  this  factor  is  dis- 
cussed in  Appendix  H.   Thus, 

"7M  *?  (25  )g  p  x  1.1 

Similar  adjustments  were  made  for  the  substitution  of 

light  water  at  508°F  for  the  heavy  water  at  room  temperature 

on  which  the  program  is  based.   Here  the  adjustment  process 

was  as  follows; 

^(moderator)     S  /?M       x    &*  H20     x  Sl±|%        t  A  D2O 

'AH20  V1+V°V*    ^DaO 


Appendix  F  lists  all  cross  sections  required  to  calculate 
both  the  cores  for  comparison  with  experiment  and  cores  for 
our  parametric  study.   Cross  section  information  was  obtained 
from  BNL-325s>  Neutron  Cross  Sections,  and  modified  for  Maxwell- 
Bolt  zmann  distribution  and  additionally  for  temperature  where 
necessary.   As  pointed  out  in  the  Introduction,  one  of  the  ad- 
vantages of  the  Fourier  program  is  its  ability  to  use  both  a 
Fermi  age  and  a  diffusion  age  and  in  this  manner  effectively 
treat  the  fast  neutrons  as  two  groups.   A  rather  complete  sur- 
vey of  the  literature  was  made  to  determine  whether  any  theoreti- 
cal work  had  been  done  on  defining  age  in  this  way  for  light 
water,  and  it  appeared  that  it  had  not.   The  most  authoritative 
experimental  results  were  computed  by  Kouts,  et.al.  at  Brookhaven 
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based  on  the  Fermi  Age  theory.   These  results  indicated  that 
a-   31  cm2  for  light  water  of  density  1  gm/cm3  for  water  to 
metal  ratios  from  one  to  three.   It  was  therefore  decided  to 
use  this  number  as  the  Fermi  age  and  to  set  the  diffusion  age 
input  equal  to  zero  thus  eliminating  this.   The  resulting  cal- 
culation performed  by  the  program  is  actually  a  one  group  con- 
tinuous slowing  down  solution.   Appendix  F  includes  a  fuller 
description  of  this  literature  survey. 

In  addition  to  the  neutron  characteristics  of  the  mate- 
rial components  of  the  reactor,  Powell's  program  also  requires 
the  dimensions  of  the  reactor  as  inputs.   Since  uniform  reflec- 
tors and  cylindrical  geometry  were  assumed  for  all  cases,  the 
only  dimensions  required  were  core  height  and  radius.   A  series 
of  set  core  heights  was  assumed;  for  given  height,  power  out- 
put, and  cell  geometry,  the  radius  is  determined  by  the  heat 
transfer  characteristics  and  temperature  limitations  of  the 
materials  and  the  overall  flux  distribution. 

A  method  of  computing  core  radius  on  the  basis  of  cell 
geometry  and  materials  for  given  power  output  and  core  height 
was  developed  from  the  functions  derived  by  Rohsenow,  Lewins, 
and  Barger  *  for  temperature  distributions  in  a  core.   This 
development, which  is  carried  through  in  detail  in  Appendix  G, 
assumes  that  the  optimum  configuration  is  that  in  which  the 
temperature  limitations  for  both  bulk  boiling  in  the  water  and 
melting  of  the  fuel  element  are  just  met  at  the  hottest  points 
in  these  media.   A  third  possible  criterion  might  have  been 
that  the  melting  point  of  the  cladding  should  Just  be  reached; 


1.  Rohsenow,  Lewins,  and  Barger,  op.clt.,  and  Appendix  G- 
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the  temperature  of  the  hottest  point  in  the  cladding  was 
checked  in  each  calculation  and  found  to  be  well  "below  the 
melting  temperature;  thus  this  was  not  a  limiting  factor. 

This  method  offers  considerable  advantages  in  a  study- 
such  as  this.   By  using  it,  it  is  possible  to  design  a  series 
of  cores  of  widely  varying  geometry  all  of  which  will  be 
capable  of  developing  the  same  amount  of  thermal  power  while 
Just  meeting  an  identical  set  of  temperature  limitations. 
This  allows  meaningful  comparisons  between  the  different  com- 
binations of  parameters. 

The  process  carried  through  on  each  of  the  cores  con- 
sidered was  as  follows? 

1.  Characteristics  were  selected  from  the  ranges  of 
parameter  variations  listed  in  the  Introduction. 
These  were:   rod  diameter,  moderator  to  fuel 
volume  ratio,  cladding  material  (stainless  steel 
or  zirconium),  and  core  height. 

2.  Using  these  characteristics,  adjusted  cross  sec- 
tions were  computed. 

3»  Using  the  same  characteristics,  a  core  radius  was 
computed;  since  this  was  the  boundary  of  a  uni- 
form lattice  of  cells  each  of  known  geometry, 
the  number  of  cells  and  hence  the  total  amount 
of  fuel  was  also  known. 

4.  Using  the  results  of  steps  (2)  and  (3),  as  input 
data,  a  Whirlwind  computation  using  Powell °s 
Fourier  program  was  made  to  find  the  critical 
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mass  of  U-235  required. 

5.  The  critical  mass  of  U-235  was  compared  to  the 
total  mass  of  uranium,  found  in  step  (3),  giving 
the  required  enrichment  of  the  uranium  for  critical- 
it  y  in  that  core. 

6.  Cores  were  compared  on  the  basis  of  amount  and 
enrichment  of  uranium  required. 

It  will  he  noted  that  enrichment  is  required  in  order 
to  find  the  disadvantage  factor,  and  therefore  in  order  to 
compute  the  adjusted  cross  sections;  it  is  also  required  to 
determine  the  amount  of  U-238  present  and  hence  the  U-238  con- 
tribution to  the  "absorber"  cross  sections.   A  solution  can 
thus  only  be  obtained  by  an  iterative  process  in  which  an  en- 
richment is  assumed  in  step  (2)  and  checked  in  step  (5)» 
However,  since  the  range  of  enrichments  Involved  is  very 
small,  from  about  one  to  three  percent,  and  the  solution 
converges  very  rapidly,  only  a  few  cores  had  to  be  run  more 
than  once,  and  none  more  than  twice. 
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III. 

RESULTS 


TABLE  I 
Total  Mass  of  Uranium  (lO^  gm) 
Stainless  Steel  and  Zirconium 
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Rod 
Diameter 

Core 

Height 

Vw/ 

Vu 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

Oil" 

8» 

6» 

,96 

4' 

1.043 

1.10 

1.19 

0.2" 

8» 

4.37 

4.57 

4.66 

4.76 

4,88 

4.92 

6» 

4.66 

4.85 

4.92 

4,96 

5,05 

ip 

4.89 

5.02 

5.13 

5.17 

5.21 

0.4" 

8« 

19-8 

20.2 

20.5 

20.7 

20.7 

6' 

20.35 

20.6 

20.9 

21.0 

21,2 

ip 

20.9 

21.1 

21.3 

21.4 

21.8 

0.6" 

8» 

46.1 

47.0 

48.0 

48.5 

48.6 

6« 

47.0 

47.6 

48.5 

49.O 

40 

48.2 

48.8 

49.3 

50.2 

TABLE  II 
Stainless  Steel  Clad 
U-235  Critical  Mass  (102*  gm) 
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Rod 

Diameter 

Core 
Height 

Vl/V0 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

0.1M 

8« 

6« 

3.72 

fyi 

6.44 

4.25 

3.60 

0.2" 

8' 

12.6 

ll.l 

8.85 

7,94 

7,82 

8.45 

6« 

11.23 

8.77 

7.61 

7.65 

8.21 

4. 

10.7 

8.76 

7. 80 

7.88 

8.44 

0.4" 

8» 

35.1 

30.3 

26.0 

24.5 

27.26 

6« 

3^.6 

29.9 

25.O 

24.1 

27.13 

40 

30.7 

28.1 

25.0 

24.2 

27.51 

0.6M 

8° 

55.4 

51.6 

50.1 

53.2 

59.9 

6' 

51.2 

50.6 

49.2 

50.8 

4« 

51.5 

49.7 

48.8 

52.22 
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TABLE  III 


Zirconium  Clad 
M25  x  10**  gn 


Rod 
Diameter 

Core 
Height 

Vl/Vo 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

0.1" 

8« 

6» 

1.88 

4» 

2.89 

2.04 

1.80 

0.2* 

8« 

7.22 

6.53 

5.60 

5.37 

5-56 

5.94 

6» 

6.64 

5.86 

5.63 

5.81 

6.20 

4. 

7.19 

6.22 

6.05 

6.14 

6.56 

o.4M 

8» 

25.4 

21.56 

18.5 

17.65 

20,37 

6» 

26.06 

22.46 

19.1 

18.67 

21.86 

41 

27.59 

23.4 

19.7 

19.38 

22,3^ 

0.6M 

8« 

43.06 

42.8 

43.0 

44.95 

51.83 

6» 

43. 81 

43.0 

43.2 

44.13 

4« 

44.81 

44.1 

44.02 

45.11 

TABLE  IV 
Enrichment 


Stainless  Steel  Clad 
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Rod 
Diameter 

Core 
Height 

Vl/Vo 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

0.1M 

8" 

6« 

3.88 

«M 

5.62 

3,75 

3.02 

0.2" 

8» 

2.88 

2.43 

1.90 

1.67 

1.60 

1.70-5 

6» 

2.41 

1.81 

1.55 

1.54 

1.625 

4,1 

2.19 

1.74 

1.52 

1.525 

1.6 

o.4M 

8< 

1.775 

1.50 

1.26 

1.182 

1.32 

6« 

1.70 

1.45 

1.20 

1.15 

1,28 

4» 

1.465 

1.33 

1.172 

1.13 

1.26 

0.6" 

8« 

1.20 

1.10 

1.045 

1.10 

1.23 

6« 

1.09 

1.06 

1.015 

1.035 

4' 

1.07 

1.02 

0.99 

1.04 

, 


TABLE  V 
Enrichment 


Zirconium  Clad 


37- 


Rod 
Diameter 

Core 
Height 

Vl/V0 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

0.2" 

8» 

1.65 

1.425 

1.20 

1.13 

1.14 

1.205 

6» 

1.43 

1.21 

1.145 

1.17 

1.225 

Jp 

1.47 

1.24 

1.18 

1.185 

1.26 

0.4" 

8» 

1.28 

1.07 

.90 

.85 

.985 

6' 

1.285 

1.09 

•  915 

.89 

1.03 

41 

1.32 

1.11 

.925 

.905 

1.035 

0.6" 

8' 

•  935 

.91 

.895 

.927 

1.065 

6« 

•  932 

•  905 

.890 

.905 

4» 

•  930 

.905 

.892 

.900 

TABLE  VI 
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Core  Diameter  to  Height  Ratio 


Rod 
Diameter 

Core 

Height 

V/vA 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

0.1" 

8« 

1 

6" 

J*i 

.60 

.66 

.754 

0.2" 

8« 

.332 

•  37 

.422 

.49 

.54 

.585 

6» 

•  577 

.68 

.765 

.84 

.917 

l±t 

1.078 

1.26 

1.41 

1.56 

1.70 

0.4" 

8» 

.705 

•  778 

.905 

1.01 

1.11 

6' 

1.11 

1.21 

1.41 

1.57 

1.74 

/p 

2.03 

2.24 

2.60 

2.91 

3.20 

0.6" 

8» 

1.10 

1.18 

1.38 

1.55 

1.70 

6' 

1.71 

1.84 

2.15 

2,41 

J}l 

3.15 

3.55 

3.95 

4.45 

TABLE  VII 


Stainless  Steel 
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Rod 
Diameter 

Core 
Height 

vi/vfe 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

0.1M 

8« 

6« 

*■ 

1.16 

1.15 

1.12 

0.2M 

8' 

1.13 

1.13 

1.11 

1.10 

1.09 

1.08 

6« 

1.13 

1.11 

1.06 

1,05 

1.G4 

fci 

1.15 

1.12 

1.08 

1.06 

1.05 

0.4M 

8« 

1.19 

1.19 

1.13 

1.11 

1,09 

6« 

1.08 

1.08 

1.10 

1,09 

1,07 

*• 

1.12 

1.11 

1.09 

1.08 

1,07 

0.6M 

8' 

1.11 

1.11 

1.09 

1,08 

1.07 

6' 

1.08 

1.11 

1.06 

1,05 

4.1 

1.0? 

1.09 

1.06 

1.05 

TABLE  VIII 
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Zirconium 


Rod 

Diameter 

Core 

Height 

▼l/Vo 

1.5 

2.0 

3.0 

4.0 

5.0 

6.0 

0.1" 

8" 

6« 

41 

1.18 

1.16 

1.15 

0.2M 

8» 

1.14 

1.15 

1.13 

1.11 

1.10 

1,09 

6« 

1.14 

1.13 

1.12 

1.11 

1,09 

J+o 

1.17 

1.17 

1.15 

1.13 

1.10 

0.4M 

8» 

1.20 

1.18 

1.10 

1.06 

1.05 

6' 

1.21 

1.19 

1.11 

1.08 

1.07 

1*1 

1.22 

1.19 

1.10 

1.08 

1.07 

0.6" 

8« 

1.12 

1.11 

1.09 

1.08 

1.07 

6" 

1.13 

1.11 

1.10 

1,07 

4.1 

1.13 

1.10 

1.09 

1.05 
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IV. 
DI30U33I0N  OF  RESULTS 

The  results  of  the  calculations  for  the  series  of  400 
megawatt  cores  described  in  the  Introduction  are  shown  in 
the  Results,  and  total  uranium  mass  and  enrichment  are  plotted 
in  Figures  2,  3  and  4.  Both  the  total  mass  of  uranium  and 
the  required  critical  enrichment  are  shown,  since  the  uranium 
cost  Is  a  function  of  both. 

Total  mass  of  uranium  is  a  function  solely  of  the  heat 
transfer  properties  of  the  combination  of  parameters  under 
consideration.   It  is  a  strong  function  of  rod  diameter  and 
a  weak  function  of  core  height  and  moderator  to  fuel  volume 
ratio.   For  a  given  core  height  the  allowable  heat  generation 
rate  in  the  central  rod  is  almost  independent  of  rod  diameter. 
This  is  due  to  the  low  thermal  conductivity  of  the  oxide 
fuel  (0.95  Btu/hr.-ft.-°F)  compared  to  those  of  stainless 
steel  (10.6)  or  zirconium  (8.15).  As  a  result,  almost  the 
entire  temperature  drop  between  the  4800°F  fuel  centerllne 
and  the  water  at  less  than  628°F  is  taken  across  the  fuel 
region.   Using  the  methods  of  Rohsenow,  Lewins,  and  Barger,— 
a  sample  temperature  distribution  was  computed  for  the 
point  of  maximum  fuel  centerllne  temperature  in  a  O.k   inch 
diameter  rod.   It  is  shown  in  Figure  1,  in  which  the  magni- 
tude of  the  drop  in  the  fuel,  in  this  case  97#  of  the  total 
drop,  can  be  seen  and  compared  with  those  in  the  clad  and 
across  the  film. 

1.  nohsenow,  Lewins,  and  Barger,  op.  cit. 
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CLAD 


WATER 


0,1  0.3L        03        OS         0,5       0*C 

RADIUS    (V)   (si*,.) 


Figure  1 

Assuming  that  the  flux  in  the  fuel  can  be  expressed  as:' 

P  -  AI0(iftr)  (1) 

then  the  heat  generation  pattern  will  be  the   same: 

q(r)   =  AI0(^0r)  (2) 

where  q(r)   is  the  volume  rate  of  heat  generation  at  radius  r. 

The  average  heat  generation  rate  for  the  whole  rod,  q,    oan  be 
expressed  as: 

q-TTRo2  =    yR°  2  7fr  AI0    (  #>r)dr 


1. 


and 


_     2JTA   .    Ro  x    (^  Ro) 
q(r)  =     ^ftp  R0  t      ,y     v 

?m%  Ro)  x°  {«or) 


(3) 


The   total  heat  generation  Q(r),    inside  any  radius,    r,   will  be: 


1.    Glasstone   and  Edlund,    9«57 


*3. 


Figure  2 
Total  Masa  of  Uranium 


X  3  H  S 

Moderator  to  Fuel  Ratio,  ^i/y 


Figure  3 

Enrichment  vs  Moderator  to  Fuel  Ratio 
for  Stainless  Steel  Clad 


kk. 


a.        3        f        5" 

Moderator  to  Fuel  Ratio,  Vi/v0 


Figure  4 

Enriohment  vs  Moderator  to  Fuel  Ratio 
for  Zirconium  Clad 
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3L        3  H  & 

Moderator  to  Fuel  Ratio,  Vi/y 
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«    ).    ■       j    2  7Tr       §Jo^o  I0(7<or)dr 

=  7rq  *o  r  i     (^0  r)  (If) 

i5T^R0) 

For  steady  state  conditions,  all  this  heat  must  be  transferred 
across  the  circumference  at  r.   Thus 


Q(r)  «  Tr,q  *? '    rl!(%r)=  -k.2im.dt/dr 

dt/dr  =       Zj^ T   •     J-  •    *l<ffor>  (5) 


and  integrating? 


V*<V  ■"*  -  ^Jfoo'fe  x  [o<*0Ho>-i| 


2 
Since  the  total  heat  generated,  Q,  is  q  TTRq  , 

dt    =    Tffl^   ^ 


The  third  factor  in  equation  (7)  is  very  nearly  constant  oyer 
the  range  of  ^QR0  involved,  ranging  from  0.499  at  %R0  =  0.05 
to  0.459  at  *JVoRo  s  0.40.   Thus,  since  2  7TK  is  constant,  for 
constant  4t,  Q,  is  the  same  regardless  of  the  value  of  Rq, 
and  very  nearly  the  same  power  per  unit  length  is  developed 
in  the  central  rod  regardless  of  its  diameter. 

Figure  5  shows  the  effect  of  increasing  the  radius  of 
the  core  for  the  same  central  rod  heat  generation,  assuming 


(6) 


5*r  l_     *L<7K>Ro)J 

A+     m       a     x  -3: x  |*o    (  ?fo  Ro)-l|  (7) 
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Figure  5 


a  J0  overall  flux  distribution.   The  average  to  maximum  heat 
generation  ratio  is  reduced  slightly,  due  to  the  large  num- 
ber of  rods  near  the  outer  circumference  which  operate  at 
low  powers.   This  difference  is,  however,  small;  the  ratios 
vary  between  about  0.5  and  0.6. 

The  two  influences  that  will  tend  to  increase  the  core 
radius  for  constant  power  output  are: 

1.  Increasing  the  volume  ratio  of  moderator  to 
fuel. 

2.  Decreasing  the  core  height. 

The  result  is  the  pattern  observed  in  Figure  2  and  in 
Table  I  of  the  Results:  since  very  nearly  the  same  power  per 
unit  length  can  be  developed  in  the  average  rod  regardless 
of  rod  diameter,  very  nearly  the  same  number  of  rods  is  re- 
quired at  a  given  core  height  regardless  of  rod  diameter, 
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and  the  mass  of  uranium  increases  roughly  as  the  square  of 
the  rod  diameter.   In  the  range  of  core  parameters  investi- 
gated, which  represents  about  the  range  of  practicality,  the 
number  of  rods  for  given  core  height  varied  through  about 
20#.   And  since  increase  of  the  moderator-fuel  ratio  increases 
the  diameter,  a  slightly  larger  number  of  rods  is  required, 
and,  therefore,  the  mass  of  uranium  is  slightly  increased. 
Similarly,  a  reduction  in  core  height  demands  an  Increase 
in  the  number  of  the  shorter  rods,  resulting  in  a  diameter 
increase,  and,  again,  a  lowering  of  the  average  to  maximum 
heat  generation  ratio  and  a  slight  increase  in  the  mass  of 
uranium. 

The  variation  of  enrichment  with  rod  diameter  and 
moderator  to  fuel  ratio  is  a  somewhat  more  complex  phenomenon. 

The  effect  of  leakage  is  very  small,  for  two  reasons* 

1.  Both  Fermi  age,^  ,  and  the  squared  diffusion 
length,  L2,  are  small  for  water  moderated  lat- 
tices.  Room  temperature  values  are  31  ©hi2  and 

4  cm2  respectively;  both  increase  with  a  decrease 
in  water  density,  so  that  at  the  test  conditions, 
the  value  s  are  ^9 • 5  cm2  and  5  cm2 . 

2.  Due  to  the  heat  transfer  requirements,  the  cores, 
of  whatever  proportion  of  length  to  diameter, 

are  relatively  large  and  the  buckling  consequent- 
ly small.   For  the  cores  considered,  B2  computed 
from  the  core  dimensions  byi-1- 

1.  (Jlasstone  and  Edlund,  7*60 


if  9. 

B2=  (^y  t  ^y 

varied  between  3  x  kH*  and  9  x  10"/* 
From  these  values  the  non-leakage  probability  can  be  com- 
puted approximately  from  the  modified  one  group  theory  dis- 
cussed In  the  Introduction,  as  the  product  of  slowing  down 
and  thermal  diffusion  non-leakages: 

slowing  down  non-leakage       ~     e"B 

thermal  diffusion  non-leakage  s  1/(1  *  L  B  ) 
total  non-leakage  a  e-B2^ /(I  *  L2B2)    (8) 

For  the  range  of  buckling s  and  the  appropriate  age  and 
diffusion  length  values,  the  non-leakage  probability  computed 
in  this  way  for  the  cores  considered  varies  between  0.95 
and  0,985.   This  is  so  small  a  range  as  to  be  nearly  constant. 

The  machine  calculation  gives  as  a  supplementary  result  the 
value  of  k00  required  for  criticality;   these  are  shown  in 
Tables  YII  and  VIII  of  the  Results;  it  will  be  seen  that  they 
vary,  for  stainless  steel  cladding,  from  about  1.06  to  about 
1.13*   Since 

k  (effective)  =  1  -  k00  x  non-leakage  (9) 

the  non-leakage  computed  by  this  somewhat  more  sophisticated 
method  varies  in  the  narrow  range  from  0.88  to  0.94^  still  very 
nearly  constant. 

This  means  that  leakage  alone  will  not  explain  the 
strong  variations  found  in  enrichment,  which  must  then  depend 
on  the  effects  of  the  cell  parameters.   These  effects  can  be 
explained  qualitatively  with  the  use  of  a  somewhat  simplified 
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model,  as  described  below. 

As  discussed  in  the  Introduction,  k00  is  the  product 
of  four  factors;  a  fifth  was  introduced  to  bring  the  computer 
results  Into  agreement  with  experiment,—1-  so  that 

k0o  =  l.l^pf 
Of  these  factors  if  is  a  constant  f or  U-235  fuel 

k00  =  1.1  x  2.08  £pf 
Of  the  remaining  three  factors,  £   and  p  were  computed  in  the 
work  preliminary  to  entering  the  machine  calculation—  and 
f  was,  In  effect,  computed  by  the  machine.   Thus, 

f=    1 — x  ^oo  =  0.i*36  —  (10) 

1.1  x  2.08    £  p  gp 

By  equation  (8), 

k00  s  1/ (total  non-leakage) 

And,  using  an  average  value  of  0.90  for  non-leakage, 


0.90     ft,  cry 


(11) 


90  -  ev  sp 

The  thermal  utilization,  f,  is  also  definied  as  the 
ratio  of  absorption  in  the  fuel  to  total  absorption  in  fuel, 
U-238,  and  moderator; 

f  s  To  ^a25    ffo     (12) 

vo   ^a25  ?o     *  V0    ^~a28  tfo  ?  ?i   Wal  fl 

Algebraic  manipulation  of  this   gives 

1  =     Vo  ^a25^o  +  Vp  ^a28  ffo     +     Vx  ^al  (fl 
f         vo  ^a25  ^0       *o  ^a25  f  o  vo   ^"a25  ^o 

1  =  1+     ?28i^     t^,xJlx^AL  (13) 

f  M25    (T25  V0      ^o         ^a25 


1.  See  Appendix  H. 

2.  See  Procedure  and  Appendix  D. 
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For  low  enrichments, 

enrichment  =  1*25/(^5  +  N28)— N25/N28 
And  for  enrichments  in  a  narrow  range  around  2%   and  cladding 
volume  much  smaller  than  water  volume, 

al/^a25  —  constant  ^  0.1 


Substituting  into  equation  (13) 


1-14  0.1  x  v] 

f "  L  ▼c 


+  ^.15  x  10'3        (ik) 
enrichment 


ro  <fca 

And  from  equation  (11) 

T=      A      =2-06"  (15) 

Equation  (15)  may  he  plotted,  as  in  Figure  6a  from  the  com- 
puted data  on  fast  fission  factor  and  resonance  escape  proba- 
bility, which  are  functions  only  of  rod  size  and  moderator 
to  fuel  ratio. 

The  two  terms  of  equation  (1*0  can  be  separated;  the 
first  term  is  a  straight  line  of  slope  0.1  with  increasing 
volume  ratio,  slightly  modified  by  the  disadvantage  factor, 

,  which  was  also  computed  in  the  preliminary  calcu- 


<Pv?, 


o 

latlons  for  each  core  and  increases  slowly  with  increasing 
rod  diameter  and  volume  ratio.   This  is  also  plotted  in 
Figure  6a.   The  difference  between  the  two  lines  so  plotted 
for  each  rod  size  must  be  the  second  term  of  equation  (1*0, 
and  therefore  proportional  to  the  reciprocal  of  enrichment. 
It  will  be  seen  from  Figure  6b,  in  which  these  differences 
are  plotted  for  the  rod  diameters  against  the  moderator  to 


Figure  6 

Effect  of  Cell  Parameters  on 
Enrichment  (Equations  (14)  and  (15)) 
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Moderator  to  Fuel  Ratio  ,  Vi/y 
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fuel  ratio,  that  the  difference  reaches  a  maximum,  represent- 
ing a  minimum  enrichment,  for  each  roa  diameter,  and  that 
these  maxima  are  reached  at  a  different  moderator  to  fuel 
ratio  for  each  rod  diameter,  which  increases  with  decreasing 
rod  diameter. 

The  same  pattern  can  be  seen  in  the  results  of  the  more 
accurate  machine  computation,  as  shown  in  Figures  3  and  k   and 
in  Tables  IT  and  V  of  the  Results.   Optimum  enrichments  are 
found  for  each  rod  diameter  at  moderator  to  fuel  ratios  which 
increase  with  a  decrease  in  rod  diameter.   If  leakage  is  as- 
sumed small  and  constant,  the  enrichment  depends  on  the  cell 
effects  shown  in  Figure  6a.  The  product  £-p  is  greater  for 
larger  rod  diameters,  and  increases  with  decreasing  slope 
as  the  moderator  to  fuel  ratio  increases.   HoweTer,  since 
water  is  a  strong  absorber  of  neutrons,  the  absorption  effect 
increases  steadily  as  more  water  is  added,  tending  to  resist 
the  advantage  of  increased  £p.   In  addition,  as  rod  size  In- 
creases, the  disadvantage  factor  effect  becomes  larger,  also 
tending  to  resist  the  effect  of  £p.  As  a  result,  an  opti- 
mum moderator  to  fuel  ratio  is  reached,  beyond  which  a  fur- 
ther increase  will  result  in  an  increased  enrichment.   Due 
to  the  earlier  leveling  off  of  £p  and  the  greater  effect  of 
disadvantage  factor,  both  of  which  can  be  clearly  seen  in 
Figure  6a,  this  optimum  occurs  at  lower  values  of  the  volume 
ratio  for  larger  rods. 

That  this  effect  is  nearly  independent  of  core  geometry 
at  this  power  is  shown  in  Figures  7  and  8a  in  which  enrichments 
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are  plotted  for  rod  diameters  and  core  heights  against  the 
diameter  to  height  ratio  of  the  core.   If  core  geometry  had 
a  significant  effect,  it  would  be  shown  by  the  existence  of 
a  minimum  enrichment  near  a  diameter  to  length  ratio  of  one. 
Although  slight  variations  exist  in  the  best  enrichment  ob- 
tainable with  a  given  rod  diameter,  they  are  significant 
mainly  in  that  they  are  so  small.   As  indicated  by  Figures 
7,  8a  and  8b,  an  optimum  enrichment  can  be  obtained  for  a 
given  rod  diameter  which  is  independent  of  the  proportions 
of  the  core  over  a  wide  range. 

Although  the  volume  of  cladding  is  small  compared  to 
those  of  fuel  and  water,  its  significance  is  seen  in  the  com- 
parison between  the  stainless  steel  and  the  zirconium  clad 
fuels,  shown  in  Figures  3»  k,  7  and  8.   Due  to  the  much  smal- 
ler neutron  absorption  in  zirconium  (  ^a  s  0.00526  compared 
*°    5*o  °  0«1535  for  stainless  steel),  critlcallty  can  be 

attained  with  considerably  lower  enrichments.   That  the 
optimum  volume  ratios  are  not  far  removed  from  those  found 
for  the  stainless  steel  clad  oores  is  due  to  the  lowering  of 

<Ea25  wi*h  "the  decrease  in  enrichment  from  about  2%   to  about 
1#.   Thus  ^"qt/  ^"a25  "the  sloPe  of  tne  lower  group  of  lines 
in  Figure  6b  before  the  disadvantage  factor  correction  is  ap- 
plied, is  not  greatly  changed,  and  the  optimum  enrichments 
occur  at  roughly  the  same  values  of  Vi/V0. 

It  should  be  observed  that  the  minimum  enrichment  does 
not  necessarily  Indicate  the  optimum  core.   The  uranium  cost 
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is  a  function  both  of  enrichment  and  the  total  quantity  of 
uranium.   Since  the  total  quantity  of  uranium  goes  up  rapidly 
with  the  rod  diameter  and  the  enrichment  is  reduced  only 
slowly.  Atomic  Energy  Commission  price  figures  released  in 
December  of  1956—'  indicate  that  the  Increased  amount  of 
uranium  required  for  the  larger  rods  will  overweight  the  low- 
er unit  cost  due  to  the  lower  enrichment,  making  the  fuel 
cost  less  for  the  0.2  inch  rods  than  for  the  0.6  inch.   This 
effect  will  be  reinforced  by  the  greater  amount  of  cladding 
material  required  to  place  an  equal  thickness  of  clad  around 
the  larger  rods,  since  the  number  of  rods  is  roughly  the  same. 
The  cost  of  rod  fabrication  and  of  the  pressure  vessel  are 
also  significant  items. 

Finally,  the  enrichments  computed  are  those  required 
to  reach  initial  criticality  only.   In  order  to  compensate 
for  the  burnup  effects  of  operation  at  power  over  a  long  period 
of  time,  these  enrichments  would  have  to  be  increased  by  a 
factor  which  depends  on  the  lattice  parameters.  A  method  of 
finding  this  factor  is  derived  in  Appendix  J.   The  cores  com- 
puted in  the  thesis  would  be  the  basis  of  such  a  calculation. 

1.  H Enriched  Uranium  Prices*,  Nucleonics  lM,    12,  R-2  (Dec.  1956) 
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V. 

CONCLUSIONS 


I.  Relating  to  methods  of  calculation. 

1.  Minimum  core  size  can  "be  computed  for  a  specified 
power  output  and  specified  unit  cell  parameters. 

2.  The  Whirlwind  program  written  by  J.  R.  Powell  using 
a  Fourier  method  of  solution  for  critical  mass  of  a 
cylindrical,  highly  enriched,  homogeneous,  heavy 
water  moderated  reactor  can  he  adapted  for  use  In  the 
cylindrical,  low  enrichment,  heterogeneous,  light 
water  case,  and  experimental  correlations  were  ob- 
tained for  the  adjusted  program. 

3*  The  factor  by  which  critical  enrichment  of  a  particu- 
lar core  must  be  Increased  to  attain  a  specified  power 
time  of  operation  can  be  computed  as  a  function  of  the 
critical  amount  of  U-235  and  the  cell  parameters. 
II.  Relating  to  the  series  of  Jj-00  megawatt  thermal  output 

cylindrical  pressurized  light  water  reactor  cores  studied. 

1.  Total  mass  of  uranium  required  is  a  function  of  the 
heat  transfer  characteristics  of  the  unit  cell;   it 
is  approximately  proportional  to  the  square  of  the 
rod  diameter. 

2.  Enrichment  required  for  initial  criticality  is  princi- 
pally a  function  of  the  unit  cell  parameters  and  is 
nearly  Independent  of  the  core  size  and  shape,  between 
diameter  to  length  of  core  ratios  of  0.8  and  3» 

3.  For  each  rod  diameter  and  cladding  material  an  optimum 
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moderator  to  fuel  volume  ratio  exists  at  which  the 
critieal  enrichment  is  a  minimum.   For  the  diameters 
tested;,  these  were  as  follows? 

Rod  Diameter  Optimum  Vi/V0 

Stainless  Steel  clad   Zirconium  clad 

0.6  inch  3„0  3.0 

0.4  inch  3„75  4.0 

0.2  inch  4.5  4.5 

4.  An  optimum  rod  diameter  at  which  minimum  critical  en- 
richment is  lowest  exists  for  rods  clad  with  each 
cladding  material.   For  zirconium  clad  rods9  this  oc- 
curs at  about  0.45  inches  rod  diameter'for  stainless 
steel  clad,  at  slightly  over  0.6  inches. 

Rod  Diameter  Minimum  critical  enrichment 

Stainless  Steel  clad   Zirconium  clad 

0.6  inch 

0„4  inch 

0.2  inch 

5«  Critical  enrichments  are  lower  for  zirconium  clad  rods 
than  for  stainless  steel  clad,  of  the  order  of  1%   en- 
richment compared  to  the  order  of  2%„ 

6.  Optimum  critical  enrichment  does  not  necessarily  indi- 
cate the  optimum  core;  other  factors  which  must  be 
considered  are  total  quantity  of  uranium,  cost  data 
for  uranium^  cladding  materials,  and  fabrication  of 
fuel  elements  and  pressure  vessels,  and  the  additional 
enrichment  necessary  to  attain  a  specified  core  life. 


.99 

.890 

1.12 

.85 

1.5 

1.13 
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APPENDIX  A 
DISADVANTAGE  FACTOR 

A.  L.  Kaplan  in  a  Massachusetts  Institute  of  Technology 
S.M.  thesis  titled  "Theoretical  Studies  of  Neutron  Flux  Dis- 
tribution Expected  in  the  MIT  Nuclear  Research  Reactor"  com- 
pares the  disadvantage  factor  calculated  by  two  group  diffu- 
sion theory  based  upon  three  reactor  models.   The  first  model 
assumes  fast  neutron  flux  constant  through  th©  cell  and  no 
net  leakage  of  thermal  neutrons  across  the  outer  boundary  of 
the  cell.   This  is  identical  to  the  development  of  Glasstone 
and  Edlund.    Model  number  two  allows  the  fast  neutron  flux 
to  vary  throughout  the  cell  as  a  function  of  the  radial  dis- 
tance from  the  center  of  the  fuel  element  and  again  assumes 
no  net  leakage  of  thermal  neutrons  across  the  outer  cell 
boundary.   The  last  model  differs  from  the  second  by  assuming 
a  net  leakage  of  thermal  neutrons  such  that  the  ratio  of 
thermal  leakage  to  fast  leakage  equals  the  ratio  for  an  equi- 
valent homogenous  cell.   He  concluded  that  the  third  model 
gave  only  three  percent  increase  in  accuracy  over  the  first 
model  and  that  this  added  accuracy  did  not  justify  the  in- 
creased complexity  in  the  calculation.   This  thesis  work  by 
Kaplan  on  heavy  water,  highly  enriched  cores  and  the  close 
agreement  between  calculated  and  experimental  disadvantage 
factors  for  light  water,  low  enrichment  cores  (Table  A-I) 
indicated  that  the  one  group  diffusion  theory  in  Glasstone 
and  Edlund  was  adequate  for  this  thesis. 
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This  development  approximates  the  reactor  core  by  divid- 
ing it  into  a  number  of  identical  unit  cells  so  that  a  square 
cross  section  is  replaced  by  a  circular  cross  section  of  equal 
area  with  radius  R^  in  the  manner  shown  in  the  following  figure 
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The  equivalent  cell  then  consists  of  a  central  fuel  rod  of 
radius  RQ  with  a  com 
clad)  as  shown  below 


radius  RQ  with  a  concentric  volume  of  moderator  (water  and 
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The  Glasstone  and  Edlund  development  assume 8  slowing 
down  density  constant  in  the  moderator  and  zero  in  the  uranium. 
With  the  boundary  conditions  of  continuity  of  neutron  flux  and 
neutron  current  density  at  the  fuel  rod-moderator  boundary 
and  no  leakage  at  the  outer  cell  boundary,  the  following  equa- 
tions for  thermal  utilization  are  developed? 

1    +  !°  5^*  *  +  <E  -  1)  (1) 

vo  2M 

p  -  Ko   Ro   Ip  QjqRq)  (2) 

ip  -   2i<gj  :  Ro2)     (io(y»Ro)KiWiRt)   g0(^'Ro)ilQr'R',n    (3) 

*       2  Ro lJ70pRi)KiQir,R,)-  K8(7r»Ho)I.  fft%  )] 

By  definition  thermal  utilization  is  expressed  by 

f  = 

Eliminating  thermal  utilization  from  equations  (1)  and  (4) 

-  P-HE-1) 


where  <$^,and  f     are  the  volume  average  thermal  fluxes  in  the 
fuel  and  moderator,  respectively,  and  their  ratio  is  the  dis- 
advantage factor.   This  appendix  contains  a  sample  calcula- 
tion for  this  factor. 

In  order  to  determine  the  accuracy  of  this  method  of 
calculation,  a  series  of  experimental  flux  traverses  for 
light  water  moderated,  low  enrichment  uranium  lattices  ob- 
tained at  Brookhaven  National  Laboratory  by  H.  J.  Kouts,  et. 
al.  was  utilized  for  comparison.   Tables  A-II,  A-III  and  A-IV 
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and  Figure  A-I  are  excerpts  from  these  reports.   Table  A-I 
compares  calculated  and  experimental  disadvantage  factors 
over  the  range  of  experimental  flux  traverses  presently  avail- 
able.  It  should  be  noted  that  this  is  an  extremely  narrow 
range;  three  enrichments,  1.027#,  1.15#  and  1.30#  for  0.6 
inch  rods  and  one  enrichment,  1.02?#,  for  0.75  inch  rods. 

Agreement  between  calculated  and  measured  disadvantage 
factors  for  the  0.6  inch  rods  are  all  less  than  6%   error. 
With  increased  diameter  the  percent  error  becomes  greater. 
This  is  explained  by  the  failure  of  the  theoretical  functions 
describing  flux  within  the  rod,  the  bessel  function  IoO^o**) 
and  flux  within  the  moderator  kl0Q{%T)+  BKt(Vir)  to  match 
the  experimental  flux  traverses  for  the  0.75  inch  rods  as 
shown  by  FigureA-  3 .   Brookhaven  National  Laboratory  found 
it  necessary  to  multiply  the  "}(    of  the  fuel  rod  by  an  average 
value  of  approximately  1.5  and  the  7(    of  the  moderator  by 
1.25  to  match  theoretical  with  experimental  flux  traverses. 
When  these  adjustments  to  7(   are  made,  the  calculated  value  is 
within  four  percent  of  the  measured  value  of  the  disadvantage 
factor  for  0.75  inch  rods. 

From  this  limited  comparison  it  appeared  that  the  one 
group  diffusion  theory  with  zero  neutron  leakage  from  the  unit 
cell  would  give  sufficiently  accurate  answers  at  very  low 
enrichments  and  at  diameters  below  0.6  inch.   However,  as 
enrichment  and  rod  diameter  increase,  the  disadvantage  factors 
become  increasingly  less  accurate.   For  the  0.75  inch  case 
the  previously  mentioned  arbitrary  factors  could  be  applied; 


6k  > 


but  until  more  experimental  work  is  done  with  enrichments 

above  J>%t  larger,  undetermined  errors  will  probably  exist  in 

the  disadvantage  factors  computed  by  this  method. 

Sample  Calculation  for  Disadvantage  Factor 

Stainless  Steel  Clad,  U02  Fuel,  H2O  Moderator  of 
.79  gm/<cm3  density 

V,/v  =  2.0  Ro  =  .508  cm 

R,     =  Ro  f^f^  ^L 

v0   =ttRo2  -B1D- 

Vc    =  (.239)(Ro+.019)  ^12i_ 

VH  0  s  ttR,2  -  (T0+Ve)  1.W* 

V,    =  VH2o  ♦  Vc  JU£2_ 

1.66  VH20/Vfl   t.^  .0115  VH20/V»  ^oxo,^ 

89  vc/v ,    -~W-  .  1535  vo/V  »  -jm*L_ 

^•/    o0225 

.040^(l=enrlch. )  0Q395 

9.65   (enricho)   .1931 

#0R0   -  .211 

Figo  IF-  1.0054- 
-  I0O36 


Figure  A-I 

Intracell  Flux  Distributions  of 
Water-Uranium  Lattices 
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2.0 


-J 
UJ 


Volume  of  Moderator  4  Clad  - 

Volume  of  Fuel 
Fuel  Rods: 0? 75  Dia. 

1.027^  u  235 

Cladding,  30  Mil  Al 


Radial  Distance  r  from  Fuel  Rod  Center  (cm) 
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TABLE  A-I 
COMPARISON  OF  CALCULATED  AND  MEASURED 
DISADVANTAGE  FACTORS 


Rod  Dlam. 

Enrichment 

ftc/Tir 

Gale.  DA 

Meas.    DA 

*  Error 

0.6 

UM" 

3 

1.32 

1.40 

5.7 

0.6 

lvTo  * 

2 

1.27 

1.30 

2.3 

0.6 

1.3a* 

1.5 

1.22 

JL.295 

5.8 

0.6 

V15  * 

3 

1.32 

1.36 

3.0 

0.6 

1.15  * 

2 

1.26 

1.33 

5.3 

0.6 

1.15  * 

1.5 

1.21 

1.27 

4.7 

0.6 

1.027* 

3 

1.29 

1.32 

2.3 

0.75 

1.027* 

2 

1.33 

1.595 

16.6 

0.75 

1.027* 

2 

1.565— 

1.595 

1.8 

0.75 

1.027* 

4 

1.5 

1.73 

13.3 

0.7*5 

1.027* 

* 

2. 

1.  fin —  ■ 

1.73 

4.0    l 

1.  Kouts  defines  Vw/Vu  as  the  rati©  of  water  to  uranium  for 
'"  0.6  inch  rods  and  as  the  ratio  of  water  plus  clad  to 
uranium  for  0.75  inch  rods. 

2^  Using  adjusted  kappa  values  of  1.5  fofor  fuel  and  1.25   #, 
for  moderator. 
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TABLE  A-II 


0.6  Inch  Diameter,  1.15^  Enrichment,  Aluminum  Clad  Uranium 
Rode  in  Light  Water  (BNL  1987) 


VVU 

fir 

& 

JLl. 

1. 

1.5 

1.449 

1.127 

1.280 

.917 

2 

1.527 

1.109 

1.250 

.888 

2 

1.514 

1.124 

1.252 

.890 

3 

1.565 

1.135 

1.290 

.842 

4 

f 

1.643 

I.129 

1.285 

.793 

i^jr  .0195 

^4    =  .0115 


1.  Dividing  the  moderator  into  cladding 
and  coolant,  thermal  utilization  is 
defined  as 


f  = 


la 


JU 


7*  *m  V«  +  fa  *w  Yw  +  *r  ^v  ti* 
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TABLE  A-III 


0.6  Inch  Diameter,  1.027#  Enrichment,  Aluminum  Clad  Uranium 
Rods  in  Light  Water  (BNL  1796) 


vw/vu 

<Pm 

<Pw 

QaL 

1. 

1 

1.079 

1.271 

1.172 

.941 

1.5 

1.105 

1.381 

1.211 

.913 

2 

1.121 

1.443 

1.246 

.886 

3 

1.107 

1.505 

1.211 

.833 

4 

1.121 

1.610 

1.252 

.781 

1. 


£M  -  .4096 

^w  s  .0195 
^irf  «  .0115 


TABLE  A-IV 


0.6  Inch  Diameter,  1.3#  Enriched,  Aluminum  Clad  Uranium 
Rods  in  Light  Water  (BNL  1783)  Values  from  graph. 
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Vw/Vu 

fa 

& 

&£ 

1 

1.35 

1.125 

1.25 

1.5 

1.45 

1.15 

1.32 

2 

1.55 

1.175 

1.35 

3 

1.75 

1.2 

1.40 
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APPENDIX  B 
METHOD  OF  ADJUSTING  HOMOGENEOUS  PROGRAM  CONSTANTS  TO  OBTAIN 
CRITICAL  MASS  OF  A  HETEROGENEOUS  REACTOR 


The  disadvantage  factor  was  utilized  to  adjust  homogen- 
eous program  constants  to  give  a  critical  mass  of  a  heterogen- 
eous reactor.   Each  value  for  a  thermal  cross  section  is  ad- 
Justed  to  have  the  same  value  as  the  effective  volume  average 
value  of  the  corresponding  cross  section  in  the  heterogeneous 
reactor  weighted  by  the  relative  flux  throughout  the  unit  cell 
in  the  following  manner. 

<  (25)  =   N(25)^(25)  Vp  f0     =  N«(25)o4(25) 
Z&{   *'  V0  f  0   ffj  ft 

where  N(25)  are  the  atoms  per  cubic  centimeter  of  fuel  and 
N0(25)  are  the  atoms  per  cubic  centimeter  of  unit  cell  volume. 

N«(25)  =   N(25)  T° 
*o  +  Vl 


erefore  /         . 


Therefore 
/ 

ft 

This  factor  was  used  to  adjust  both  U-235  and  U-238  absorption 
and  transport  cross  sections. 

Light  water  cross  sections  were  adjusted  in  an  analogous 
manner 

^    (w)   =     W(w)<ri    (w)   Yi^         =  N'(w)(7f  (w) 

^a  Top0tTif, 

N°(w)   =     N(w)   YX 
vo  +  Yx 

at*  (w)  =  (y°  "LYl \  crzM 

This  factor  was  used  to  adjust  water  absorption  and  transport 
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cross  sections. 

Cladding  cross  sections  were  not  adjusted  because  the 
flux  In  the  clad  Is  very  nearly  equal  to  the  average  flux 
In  the  unit  cell.   This  Is  apparent  by  Inspection  of  Figure 
A-I.   This  fact  was  checked  by  calculation  by  comparing 
values  of  y<S)M6  with  experimental  values  of    /f*l 
where 


?«*        v°&  +  vi 


Mf* 

The  error  in  all  cases  indicates  that  by  not  adjusting  the 
cross  sections  for  cladding  a  greater  cladding  absorption 
was  being  assumed.   This  was  a  result  of  the  fact  that  the 
average  flux  in  a  unit  cell  was  higher  than  the  actual 
neutron  flux  in  the  clad.  However,  the  error  introduced  is 
necessarily  small  because  of  the  small  percentage  of  volume 
that  the  clad  occupies  in  a  unit  cell. 
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APPENDIX  C 
ADJUSTMENTS  TO  NUCLEAR  CONSTANTS 
FOR  POWELL  FOURIER  PROGRAM 
The  Powell  program  was  written  for  a  highly  enriched, 
homogeneous,  heavy  water  moderated  reactor.   The  adjustment 
of  the  cross  sections  to  obtain  a  heterogeneous  critical  mass 
is  discussed  in  Appendix  B.   Table  E-=I  lists  the  inputs  on 
the  data  tape  for  Powell's  cylindrical  Fourier  program. 

In  this  thesis  it  was  planned  to  utilize  the  program 
for  low  enrichment,  light  water  power  reactors.   However,  there 
was  no  provision  for  the  handling  of  uranium-238  in  the  pro- 
gram.  The  best  method  appeared  to  be  to  continue  t©  define 
V   as  V  &L  and  to  treat  the  uranium-238  as  an  addi- 

tional  cladding.   By  this  method  the  absorption  in  the 
uranium-238  would  appear  in  the  four  factor  formula  in  the 
thermal  utilization. 

Jb  -  v?/ 

f      s  thermal  absorption  in  U-235 

total  thermal  absorption  in  core 

The  fast  fission  factor, £  ,  and  the  resonance  escape  proba- 
bility, p,  were  combined  with  the  7[   for  uranium-235.   This 
product  gi&y  with  an  additional  arbitrary  factor  of  1.1  to 
bring  experimental  and  program  critical  masses  into  agreement, 
constituted  the  input  for  "^25"  ln  the  Powe11  program. 

The  program,  as  written,  included  density  and  atomic 
weight  for  aluminum  and  heavy  water  as  the  cladding  and  modera- 
tor, respectively.   In  order  to  use  the  program  for  other  materials. 
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it  was  necessary  to  either  modify  inputs  to  account  for  these 
constants  or  to  modify  the  program  by  changing  the  inputs  of 
the  four  drum  addresses  concerned.   Both  methods  were  utilized. 

Method  A  combined  the  effect  of  uranium-238  and  clad 
material  in  the  following  manner 

Tklv*  *4*Ke  <3) 

iV<,Wc){&-CG.S*IO~*1!i) 
%(PJ^  l-£±L-\  o£  (H,D)  (5) 


The  drum  addresses  were  modified  to  read 

DA  420/  420/  /°CkO    s  A** 

DA  42V  424/  A^p   =   t  18. 

DA  418/  418/  Af   =   *  1. 

DA  422/  422/  $tf   S   +1. 

This  method  used  the  macroscopic  cross  sections  divided 
by  Avogadro°s  number  and  modified  the  density  and  the  atomic 
number  of  the  clad  to  equal  one  so  that  when  the  computed 
N(al)  and   V^(al)M  were  multiplied  together  an  equivalent 
clad  macroscopic  cross  section  was  obtained.   Transport  cross 
sections  were  computed  identically. 

Method  B  is  a  mathematically  equivalent  method  and  was 
devised  to  circumvent  the  modifying  of  drum  addresses.  This 
was  done  in  the  following  manners 


7*. 


(   V.  +  V,      \ 


Htx0      =   teunio  />   »    (Lo^ho 


(6) 


(7) 


(8) 


Method  B  is  merely  a  variation  of  Method  A  allowing  the 
use  of  the  program  without  modification  to  existing  constants 
incorporated  in  the  writing  of  the  program.  Method  A,  when 
used  in  the  program,  gave  inconsistent  low  critical  masses 
which  led  the  authors  to  believe  that  the  drum  addresses  were 
not  being  properly  modified,.   This  was  confirmed  by  the  sub- 
sidiary outputs  of  the  program,,  The  following  page  is  a  sample 
calculation  for  UO2  fuel  and  light  water  moderator  of  density 
0.79  gm/cm3  employing  method  B. 

The  quantity  MxM  in  the  program  was  changed  from  the 
ratio  of  the  volumes  of  aluminum  to  heavy  water  to  the  ratio 
of  the  sum  of  the  uranium-238  and  cladding  material  to  light 
water.   wxM  is  used  in  the  program  to  homogenize  the  macroscopic 
cross  section  of  cladding  and  moderator. 


£(*)(— JU ,)   4.  V/arf/  V™  t  /4 )     do 
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Sample  Calculations  for 
Nuclear  Data  U02.  Fuel,  H2O  Moderator,  Pw  =  .79  (Method  B) 


28 

.04o4 

.212 

55 

.1535 

.890 

vo 

-   .810 

vo 

~   .126 

Vq  +  Vc 

=   .936 

y0  +  Vc  x. 06025  s   q0564   (A) 
I 


.3^  x  (C 
50.0  x  (C 


(.98)(e8l)f-  .126 
lo494 


/>„-. 


V0|£_   =.783 


:b; 


:o 


Afi" 


•<£a|" 
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APPENDIX  D 
COMPUTATION  FOR  FAST  FISSION  FACTOR 
AND  RESONANCE  ESCAPE  PROBABILITY 

In  order  to  account  for  the  fisaions  In  uranium-235 
and  uranlum-238  produced  by  fast  neutrons  before  slowing 
down  to  thermal  energies,  a  fast  fission  factor  is  utilized. 
It  is  defined  as  the  rati®  of  the  total  number  of  fast  neu- 
trons produced  by  fissions  due  to  neutrons  of  all  energies 
to  the  number  resulting  from  thermal  neutron  fissions.   Com- 
puted fast  fission  factors  were  calculated  by  using  Hellen°s 
empirical  formulations 


1  "*"  1*  0.875 /»w  £l  +  .288  ST  (1) 

f     Ml  VU 


Pyt   a  density  of  water 
Vw  s  volume  of  water 
?u  s  volume  of  uranium 
Vs  =  volume  of  structural  materials 
This  factor  was  checked  for  eight  critical  experimental 
cases  performed  at  the  Bettis  plant  of  the  Westlnghouse  Corp. 
for  which  a  summary  is  given  in  Table  D-I.   Errors  of  computed 
values  are  In  all  cases  less  than  1  percent  of  measured  values. 

That  fraction  of  source  neutrons  which  escape  capture 
while  being  slowed  down  to  thermal  energy  Is  represented  by 
the  resonance  escape  probability;,  p.   The  method  outlined  in 
Glasstone  and  Edlund  for  heterogeneous  reactors  was  used  with 
the  empirical  quantities  from  The  Reactor  Handbook,  Volume  I, 
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to  compute  resonance  escape  probability. 

This  method  assumes  that  neutrons  within  the  energy  range 
where  resonance  absorption  by  fuel  is  important  can  be  repre- 
sented as  a  single  energy  group,  using  effective  slowing  down 
and  absorption  cross  sections. 

A  resonance  utilization  analogous  to  the  thermal  utiliza- 
tion is  defined  as? 


f   a  rate  or  absorption  of  resonance  neutrons    /£) 
r     rate  of  production  of  resonance  neutrons 


V 


($*)r  £*0 


v0  (JlfcJ^  *?,<!,  )r£,  (3) 

Equation  (3)  assumes  that  the  leakage  in  the  resonance  group 
is  small  compared  to  absorption  and  slowing  down.   Resonance 
utilization  may  be  expressed  in  terms  of  lattice  parameters 
ass 

F  +  (E-l)  W 


where  F  and  E  for  cylindrical  geometry  are  defined  in  Appendix 
A. 

Resonance  escape  probability  is  written  in  terms  of 
i  *»  a  s  o 

p  -  exp  (-  j-i-f  )  (5) 

It  has  been  assumed  that  the  resonance  neutrons  can  be 
approximated  by  one  energy  group.  They  exist  over  an  energy 
range  greater  than  would  be  predicted  by  one  group  diffusion 
theory  because  of  the  selective  absorption  of  neutrons  as  they 
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pass  through  uranium-238  resonances  enroute  to  the  interior 
of  the  fuel  rod. 

To  determine  this  effect  average  slowing  down  and  ab- 
sorption cross  sections  are  computed  based  upon  empirical 
relations  derived  from  experiment.   The  microscopic  average 
resonance  cross  section  is  expressed  ass 

In  E2/E1 

M  =  mass  of  fuel  lump,  grams 

S  -  surface  area  of  lumpp  cm2 

A  =  empirical  constant ,  barns 

ft  =  empirical  constant,  gm/cm2 

E^  =  lower  energy  limit  of  resonance 
E2  =  upper  energy  limit  of  resonance 

Table  D-II  lists  empirical  constants  from  The  Reactor  Hand- 
book, Volume  I,  for  natural  uranium. 
For  cylindrical  rods? 

S  =  27rRph   «.  2_ 
M    /3  7TR0*h    *R0 


NA  [1+4-  4j      =      NA  Q*  jhkj 

In  E2/E1  In  E2/E1 


The  effective  resonance  slowing  down  cross  section  for 
the  moderator  may  be  expressed  bys 


Table  D-IIX  lists  the  resonance  constants  for  moderator  from 


79- 


The  Reactor  Handbook,  Volume  I. 

A  sample  calculation  sheet  for  the  calculation  of  reson- 
ance escape  probability  for  UO2  fuel  is  shown  on  the  follow- 
ing page. 

Eight  experimental  oases  were  calculated  and  comparison 
of  measured  and  computed  values  are  shown  in  Table  D-I.   The 
maximum  percent  error  from  the  measured  p  is  2.k%. 
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Sample  Calculation  for 
Resonance  Escape  Probability U O2  Fuel,  H2O  Moderator  (.79  gm/cn?  ) 

Core  No.  39 

Fuel  Rod  Radius,  R0  =  .508  cm 
Unit  Cell  Radius, R 1  =  .880  cm 

tTr?33    *  ,0353     =    -0620     *** 

.42  x  R0  =  .213     (%Ho) 

.649*  R0  s  -329     (7f,Ro) 
.649*  R,  =  .570     {%(,) 

From  Appendix  A,        F  =  1.005 

B  =  1.026 

Iw  x  .1&0   F  .  1.494  x  .140  x  1.005  =  4.22   (A) 
V0    g  .810  x  .0626  — 

fr  "      Cf)*    E J  "        =    &.22  -f  1.026J  ~1  r    .190 

rrr        T^IO  _(B) 

P  =    e  =  .790 

Fast  Fission  Factor 

1  +  .875/.w^       .288% 

Oiii -  i.o4i 

,   «   „w        x  1.494  a.  (.288)  (.126) 


=  1  ♦ 


l.lW-  (l.l)(1.04l)(.790)(2.08)   =  1.885 


81. 


TABLE  D-I 
COMPARISON  OF  MEASURED  AND  CALCULATED  FAST  FISSION 
FACTORS  AND  RESONANCE  ESCAPE  PROBABILITIES 


Fuel 
U 

Rod 
Diam. 
in. 

Clad 
Al 

Enrich. 
% 

Vw/Vu 

Calc. 
t 

1.043 

Meas. 

e 

terror 

Calc. 
P 

.859 

Meas. 
P 

terror 
P 

.6 

1.3 

3/1 

1.04? 

.4 

.850 

1.1 

U 

.6 

Al 

1.3 

2/1 

1.056 

1.061 

.5 

.796 

.79^ 

•  3 

U 

.6 

Al 

1.3 

1.5/1 

1.066 

1.072 

.5 

.737 

.720 

2.3 

U 

.6 

Al 

10.5 

3/1 

1.043 

1.050 

•  7 

.859 

.839 

2.4 

U 

.6 

Al 

1.15 

2/1 

1.056 

1.065 

.9 

.796 

.792 

.5 

UO2 

,6 

Al 

1-3 

3.05/1 

1.041 

1.039 

.2 

.816 

.810 

•  7 

UO2 

.6 

Al 

1.3 

4/1 

1.033 

1.033 

0 

.834 

.835 

2.2 

UO2 

.6 

Al 

1.3 

5/1 

1.031 

I.029 

.2 

.879 

.870 

1.0    | 
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TABLE  D-II 

1. 

EMPIRICAL  RESONANCE  CONSTANTS  FOR  FUEL — 


Material 

A  (barns) 

//  (gm/cm2) 

In  E2/Ei 

ft>m-l) 

U  metal 

9.25 

8/3 

5.6 

0.^2 

UO2 

11.51 

1.92 

7.3 

O.Jf2 

1.  Reactor  Handbook,  Vol.  I 


TABLE  D-II I 
EMPIRICAL  RESONANCE  CONSTANTS  FOR  MODERATOR— 


Material 

Density 
(gm/om2) 

barns/atom 

Vr  2. 

cn2/gm  — 

H20 

1.0 

38.5 

•583 

1.  Reactor  Handbook,  Vol.  I. 

2^  ^fe  is  multiplied  by  .88  for  U02  fuel 


r 
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APPENDIX  E 


POWELL  CYLINDRICAL  FOURIER  PROGRAM  INPUTS  AND  OUTPUTS 

Figure  E~I  defines  the  dimensions  for  the  six  radial 
and  five  longitudinal  regions  for  Powell's  program.   In  this 
thesis  two  radial  regions;  a  central  core  and  radial  reflec- 
tor and  three  longitudinal  regions;  a  central  core  and  a 
top  and  "bottom  reflector  were  used.   Thus  R]_  and  R2  were 
effectively  set  equal  to  zero,  i.e.  1  x  10~10cm.  R3  was  the 
core  radius  computed  from  heat  transfer  considerations.   The 
thickness  of  the  two  outer  reflectors  was  set  equal  to  zero 
by  letting  R3  =  R4  =  R5.   R6>  the  outer  radius  of  the  reactor, 
was  then  simply  the  sum  of  the  core  radius  and  reflector 
thickness.   In  every  case  it  was  assumed  that  a  reflector 
thickness  of  eighteen  centimeters  was  equivalent  to  an  in- 
finite reflector.   The  radial  reflector  thickness  was  set 
equal  to  eighteen  centimeters.   The  thickness  of  the  top  and 
bottom  reflectors  (number  two)  was  set  equal  to  zero.   The  top 
and  bottom  reflector  (number  one)  thickness  was  set  equal  to 
eighteen  centimeters. 

For  the  radial  region  chosen  as  the  oore,  in  this  case 
region  three;  Fermi  age,  thermal  diffusion  coefficient,  in- 
finite multiplication  factor,  and  thermal  absorption  cross 
sections  are  computed.   Thus,  the  inputs  to  these  regions 
are  normally  chosen  as  zero  except  for  the  thermal  diffusion 
coefficient  which  appears  in  the  denominator  and  must  not  have  a 
zero  value.   This  same  situation  exists  for  these  quantities 
in  the  longitudinal  properties.   The  core  properties  are 
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replaced  by  computed  numbers  as  the  program  runs  so  that  any 
number  used  to  describe  the  core  region  will  be  erased. 

The  inputs  in  the  nuclear  constant  section  are  discussed 
in  Appendix  C 

The  final  inputs  are  two  k00  guesses  which  are  used  to 
compute  values  of  the  critical  equation.   In  general,  neither 
of  these  guesses  will  satisfy  the  critical  conditions,  so 
the  two  assumed  values  of  kQo  corresponding  to  the  tw©  com- 
puted values  of  the  critical  determinant  are  extrapolated 
linearly  to  the  value  of  k0o  a"t  which  the  value  of  the  deter- 
minant is  zero.   This  process  is  continued  always  throwing 
away  the  guess  that  is  most  in  error  until  the  difference  be- 
tween M(25)n  -3  and  M(25)n  -2  Is  less  than  one  tenth  of  a  gram. 
The  guesses  need  net  be  close  to  the  answer  except  that  the 
closer  they  are,  the  less  time  will  be  required  for  convergence 

Of  particular  interest  when  making  a  parametric  study 
such  as  this  is  the  fact  that  once  the  main  program  is  fed  into 
the  computer,  a  series  of  data  tapes  may  be  run  without  read- 
ing in  the  main  program  again.   The  data  tapes  need  not  be 
complete,  but  need  only  be  cut  from  the  beginning  through  the 
numbers  which  were  changed  from  the  preceding  data  tape.   All 
constants  read  in  on  the  first  parameter  remain  in  storage 
until  replaced  by  a  new  number,  and  all  computed  values  in 
storage  are  erased  prior  to  starting  a  new  run.   Thus  when 
systematically  varying  one  constant  a  great  amount  of  time 
can  be  saved. 

Table  E-I  is  a  sample  data  tape  preparation  form  from 


85. 


the  J.  Pewell  thesis  using  data  from  run  39.   Table  E-II  is 
a  sample  of  the  form  found  convenient  for  typing  this  data 
on  the  Whirlwind  data  tape. 


Figure  E-I 

Powell  Fourier  Program 
Core  Configuration 
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TABLE  E-I 
DATA  TAPE  PREPARATION  FORM 
FOURIER  METHOD  FOR  CYLINDER 


fc  TAPE  40^-307-39 

NOT  PA 

///Dimensions 
DA  32/32/ 


LEE  AND  GARDNER 


+  0. 


4-1.0  x  10"10  Rx 
4-1.0  x  10~10  R2 


±2!hi 


±2!L1 


±2iii 


+112.5 


+243. 


■H8. 


4-0. 


4-18. 


4-  0. 


R3   (core  radius) 
R4 

R6  (outer  radius) 

Core  thickness 

Thickness  top  reflector  1 

Thickness  top  reflector  2 

Thickness  bottom  reflector  1 

Thickness  "bottom  reflector  2 


///  Radial  properties 

///  Thermal  absorption  cross  sections 

+  0. ^al 

^a2 


40. 
4-0. 


+-0. 


-  ^a3 


88 


+  .01145  £&6 
///  koo  4-  0.  k*,    1 

+  0-  k^    3 

+  0.  ]c^    4 

40.  k«o    5 

4  0.  k^    6 

///  Thermal  diffusion  coefficients   (non  zero  values) 

£i »2 


■f  1.  Do 


+  1. D^ 

■J-  1.  D^ 


///  Ages 


t-jl d6 

f  0.  r  2 

449-5  r6 
///  Nuclear  constants 

+  421.  C^a    (25)    (barns) 

+  350.  ^  f   (25)    (barns) 

+  1.885  ^     (25) 

•^2.46  V      (25) 
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3 


+   .910      ^-a  (Al)  (bams) 

£   .3^7      cr^SL   (D20)  (barns) 

•f  .616        x  =  Vol.  Al/  Vol.  D2O 

f  1. __  Fermi  fraction  of  age 

•hO. Diffusion  frac.  of  age 

f  0.   


///  Longitudinal  Properties 

■f-  .  01145  ^a  Bottom  reflector  2 
+  .01145  ^a  Bottom  reflector  1 
+  0. 


///* 


00 


■f  .01145  ^"a  Top  reflector  1 

jh  .01145  ^"a  Top  reflector  2 

4-  0.  k^  b2 

+  0.  k^  bl 


**>  tl 


■f  0.         kpO  t2 
///  Thermal  diffusion  coefficients 

t-tS D  b2 

•f  .2 D  bl 


///  Ages 


«f  1. 
t'2 -  Dt2 


+»9.5  7b2 

M9.5  T^bl 

*^9.5  #   tl 

+  *»9.5  7"  t2 


///  Initial  k^o  guesses 

JL. 
1. 


1.1    ,     k  _  1 


1.2         k^  2 


0. 


0. 
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Notes   All  numbers  must  have  -f  or  -  and  decimal  points. 
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TABLE  E-II 
SAMPLE  DATA  TAPE  PRINT  OUT  -  RUK  39 
fe  TAPE  404-306-39      LEE  AND  GARDNER 
(24,6) 


NOT  PA 

DA  32/32/ +0. 

4  1.0   X  10~10 

tl.O 

x  10-10 

49^.5 

49^-5 

4  9^.5 

4  112.5 

4-243. 

+  18. 

+  0. 

418. 

40. 

4  0. 

4  0. 

4   0. 

40. 

40. 

-h.  01145 

4-  0. 

+  0. 

*o. 

+  0. 

4  0. 

4*0, 

4  1. 

4  1. 

4  1. 

41. 

4>1. 

4-  .2 

4  0. 

+  0. 

4  0. 

4-0. 

+  0. 

449,5 

♦  421.      4 

•350. 

4  1.885 

42.46 

4  .910 

4.3V? 

4    .616 

449.5 

4  0. 

4  0. 

40. 

•450.5 

4  4.97 

♦  1.0 

4>  0. 

4  0. 

4  .01345     +  ,01145    +.01145 

4.01145        4,01145 

+  0. 

+  0. 

4-0. 

4-0. 

f-0. 

*  .2 

+  .2 

4  .2 

4>.2 

4   .2 

+^9-5 

+^9.5 

+  ^9,5 

4-49.5 

4  1. 

41. 

41.1 

41.2 

40. 

4*0. 

1  START  AT  186 
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OUTPUTS  OF  POWELL  CYLINDRICAL  FOURIER  PROGRAM 

The  program  used  in  this  thesis  eomputed  critical  mass 
and  infinite  multiplication  factor.   However,  the  printout 
format  of,  the  program  includes  several  other  numbers  which 
can  be  very  useful  in  the  evaluation  of  the  answer  received 
and  also  very  confusing  to  the  uninitiated. 

An  actual  printout  received  for  run  39  is  shown  below: 
TAPE  40^-306-=39      LEE  AND  GARDNER   05*H\2     4~l6»57 
DECIMAL 

-K  13^6l303/-00  -K17050170/-00 

+•.  304811 65  /+0  6  +.30373265/+06 

-.15789275/-01-  .72544313A02 

*. 11931752/  01  +.11916215/+01 

+  .30481165/  06  =  .15789275/=01 

+  .  3037327V  06  -.  70 792^10/+  03 

+  .16155871/-27+.16155871/-27  +ol6l5587l/-27  +.681 7410  7/+07 

+-.2900473V-00  +.49500001/+02  +.49500001/+02  +.49500001/+02 

+.l6l5587l/=27  +.16155871/-27 
+.23188179/-01  +.20502647/-01 

Interpretation  would  be  as  follows? 

I*  the  tape  number  the  404  is  the  problem  number,  the  306  the 
programmer  s  number,  and  the  39  the  run  number.   The  numbers 
following  the  name  are  the  time  and  date  of  the  run.   In  each 
number  printed,  the  numbers  after  the  slant  (/)  locate  the 
decimal  point  of  the  preceding  number  (+.119l6215/f01  =  1.1916215) 
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The  number  +. 161558 71/-2 7  signifies  zero.   The  first  eight 
numbers  are  the  following? 

n  =  2 
Time  started  Not  used 

Critical  mass-M(25)  (grams)        304811 . 65 
Value  of  critical  determinant-^   -.015789275 
Infinite  multiplication  factor-^  1.1931752 

The  next  four  numbers  represent  the  approximation  before 
the  final  number  printed  out  above  for  the  n  -  2  and  the  n  -  3 
terms  in  the  truncated  Fourier  series . 
n  =  2  (M(25))  n  -  2  (A) 

i-1  approx.  -  304811.65  gm.       i-1  approx.  =  -.015789275 


n  =  3 
Not  used 
303732.65 

-.00725^313 
1.1916215 


n  =  3  (M(25)> 

i-1  approx.  s  303732o74  gm. 


a  =  3  (A) 

i-1  approx.  =  -.0007079241 


The  following  ten  numbers  are  read  across  the  page  from 
left  to  right  and  signify  the  following? 
The  first  three  always  remain  .16155871  x  lO"2?  =  0 

(4)  core  volume  —  6817410.7  cm3 

(5)  core  thermal  diffusion  coefficient  —  .29004734  cm. 

(6)  core  age  —  49.500001  cm2 

(7)  average  reactor  age  —  49.500001  cm2 

(8)  average  Fermi  age  —  49. 500001  cm2 

(9)  average  diffusion  age  ~  .16174107  x  10~27  ■  q 

(10)  This  number  always  remains  zero  signified  by  .16174107  x  10~2? 

The  last  two  terms  are  the  number  of  atoms  per  cubic  centi- 

=24 
meter  x  10    of  clad  and  moderator,  respectively. 
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(1)  clad  —  .023188179  atoms/cm3  x   i(r2V 

(2)  moderator  —  .0205026^7  atoms/cm3  x  ltr2^ 

With  these  outputs  It  is  possible  to  compare  each  of 
the  outputs  with  inputs  to  determine  whether  the  program  is 
operating  satisfactorily.   In  addition  by  comparison  of  the 
i=l  and  i  approximations  for  each  of  the  terms  in  the  trun- 
cated Fourier  series,,  it  is  possible  to  establish  how  fast 
the  program  is  converging. 
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APPENDIX  F 
NUCLEAR  DATA 

A  survey  of  the  literature  was  undertaken  to  determine 
best  values  of  the  nuclear  data  that  would  be  required.   In 
general  the  best  sources  of  information  were  found  t©  be  in 
the  various  Brookhaven  National  Laboratory  reports  and  in  the 
Geneva  papers. 
A.  Cross  sections 

All  cross  sections  were  obtained  from  BNL=325  except 
where  noted. 
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2. 


TABLE  F 

=1 

CROSj 

(barns ) 
.330 

3  SECTIONS 

(2200  m/ 

sec) 

(barns) 

Nuclide 

(barns) 
•  293 

(barns) 

(barns) 

-24 
NxlO 

atoms/cm-^ 

H 

38.0 

.0334 

H20 

.660 

.586 

.0334 

U-235 

687 

599- 

10. 

580. 

503. 

.0474 

U-238 

2.75 

2.44 

8.3' 

.0474 

Al 

o230 

.204 

1.4 

.06025 

Fe 

2,53 

2.24 

11. 

.0622 

Ni 

4.6 

4.07 

17*5 

.00727 

Cr 

2.9 

2.57 

3.0 

.0149 

Zr 

.19 

0I68 

8.0 

* 

.0423 

1.  d^o   signifies  effective  cross  section-corrected  for  Maxwell 
"  Boltzmann  distribution. 


2.  Type  Nr.  304  stainless  steel  was  used. 
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CRO 

TABLE 

F-II 

F) 

S3  SECTIONS  (508° 

Nuclide 

ii. 

2) 

-24 
NxlO  ** 

atoms/cm3 

.0474 
.0224 

barns 

1.80 
1.80 

cm"-1- 

.0839 
.0404 

OX**- 

barns 

9.45 
9o45 

cm""1 
.440 
.220 

barns 

U-238 
U-238(UO 

U-235 

,0474 

431. 

20.4 

96.5 

4.57 

358. 

Fe 

.0622 

1.66 

.1033 

II.  12 

.694 

Ni 

.00727 

3.01 

.0219 

20.0 

.146 

Cr 

.0149 

1.90 

.0283 

3.34 

.050 

Stain  S. 

.1535 

.890 

Al 

.06025 

.151 

.0091 

1.365 

.0823 

4. 
H20— 

.0264 

.432 

. 01145 

62.9^- 

-L66 

Zr  *L 

.0423 

.1245 

0 00526 

7^95 

•  336 

655  c^ 


c?x 


2^  o^   =^(i-a) (1-4/5^);  As  3A 

2^  Experimental-Hughes,  Pile  Neutron  Research,  page  223 

4^  Density  of  water  »  0.79  gm/cm3  (508QF) 

^  Density  of  U02   s  10.07  gm/em3 

6.  Density  of   Zr  =6.4  gm/cm3 
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B.  Fermi  Age 

Early  measurements  of  Fermi  Age  were  carried  out  at  Oak 
Ridge  in  19*j4-45 ,  for  water  to  metal  ratios  in  the  range  from 
1  to  3.5o  using  rod  diameters  of  1.18,  1.10,  and  .787  inches 
and  enrichments  of  the  order  of  1%.     A  fission  neutron  source 
was  used,  and  the  rods  were  wrapped  in  cadmium  to  prevent  fur= 
ther  fissions  in  the  lattice.   The  activities  of  cadmium  cov- 
ered indium  foils  spaced  through  the  lattice  were  measured, 
and  age  to  indium  resonance  computed  by 


^-\-n 


A  correction  term  for  the  difference  between  indium  resonance 
and  thermal  was  added,  but  this  is  small  (about  1  cm^). 

These  measurements  gave  a  set  of  characteristic  curves 
showing  an  age  of  about  40  cm  at  a  volume  ratio  of  3  to  1, 
rising  sharply  as  the  volume  ratio  was  reduced  below  1.5»   Ai 
increase  also  occurred  with  reduction  of  the  rod  diameter. 

Kouts,  in  Geneva  Paper  No.  600,  "Exponential  Experiments 
in  Light  Water",  has  reassessed  these  data  In  terms  of  the  pos= 
sibility  of  distortion  through  fast  and  eplthermal  fissions. 
These  would  not  be  prevented  by  the  cadmium  wrappings  and 
would  tend  to  occur  more  strongly  at  the  lower  volume  ratios. 
They  would  have  the  effect  of  increasing  the  apparent  age. 
Kouts  states  that,  in  fact,  the  fraction  of  fissions  produced 
by  neutrons  above  cadmium  resonances  (0.^  ev)  at  a  volume  ratio 
of  3  is  about  .11  while  at  a  volume  ratio  of  about  1.5  it 
rises  to  .58;   this  includes  both  fast  fission  in  U-238  and  . 
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epithermal  in  U-235. 

o 

Kouts  made  measurements  of  M  ,  the  migration  area,  in 

low  enrichment  water  lattices  with  rolume  ratios  of  1  to  3* 
The  age  computations  which  he  felt  were  the  most  accurate 
were  based  on  the  relationship 

His  method  was  that  of  poisoning  the  lattices  with  increasing 
amounts  of  boric  acid,  computing  B2  by  differentiation  of  the 
thermal  flux  (found  by  measuring  indium  foil  activities),  com- 
puting f,  and  solving  for  M2,  considering  the  other  factors 
to  remain  constant., 

His  results  were  markedly  different  from  the  Oak  Ridge 
curves.      They  show  a  curve  very  slowly  rising    with  decreasing 
volume  ratios,   and,    indeed,    in  one  case  flat  within  the   experi- 
mental error.     Measurements  made  with  0.600  inch  rods  are  con- 
sistently greater  than  those  with  0.750  inch  rods,   but  can  still 

p 

be  considered  almost  flat,  with  a  mean  value  of  about  30  cm  . 

Ghernick  discusses  in  his  paper,  read  at  Geneva  (No.  603, 
"The  Theory  of  Uranium  Water  Lattices")  methods  of  calculating 
age,  and,  while  he  does  not  show  actual  calculations  or  tabu- 
late results,  he  does  give  a  simplified  formula 

r-   \  [1+  3,+**-] 

which  does  not  include  volume  ratio  as  a  variable  and  states 
that  "Calculated  values  of  the  age  of  fission  neutrons  to 
indium  resonance  in  water  to  uranium  mixtures  show  little 
variation  if  the  uranium  to  water  volume  ratio  is  less  than 
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unity. • . ,Fbr  equal  volumes  of  uranium  and  water,  calculated 
ages  have  run  about  20#  higher  than  the  value  for  pure  water." 

Cherniek  gives  as  the  best  current  value  for  age  to 
indium  resonance  in  pure  water,  25  cm  . 

1.2  x  25  =  30  cm2 

Also  Included  in  Oherniok's  paper  is  a  formula  attri- 
buted to  Ehrlich  and  Deutsch  of  KAPL 


(Afvx)* 


where  v  =  volume  ratio,  uranium  to  water 

~   density  of  water 
x  *  O.38  for  stainless  steel 
0.^0  for  Zr. 


r. 


7  ~-   25  (1^1)2    =  52.5  cm^ 


(1  Y-.38)2 

Deutsch  in  Nucleonics  (January  1957)  gives  a  method  for 
computing  age  in  mixtures  of  waters  and  various  metals  such 
as  Al,  Fe,  and  zr,  for  high  enrichment  reactors.   These  ages 
will  be  in  the  range  of  *f0  to  50  or  more,  but  they  are  appar- 
ently based  on  significant  volumes  of  cladding  material  and 
insignificant  volumes  of  uranium.   These,  then,  do  not  fit  our 
problem.   Because  of  the  similarity  between  this  method  and 
the  formula  ascribed  to  Deutsch  by  Cherniek  it  seems  possible 
that  that  formula  might  well  have  been  intended  for  use  with 
high  enrichments  also,  and  the  "vM  in  it  might  have  been  in- 
tended to  express  total  metal  to  moderator  ratio,  being  largely 
made  up  of  cladding  volume. 


101. 


Finally,  J.  E.  Wilkins  et.  al.,  In  Geneva  Paper  No.  59?', 
"Status  of  Experimental  and  Theoretical  Information  on  Neutron 
Slowing  Down  Distributions  in  Hydrogenous  Media",  discuss  com- 
putation methods  for  strongly  absorbing  media  in  which  thermal 
diffusion  area  is  small  compared  to  age,  and  compares  results 
to  Kouts °  data.   The  calculated  line  follows  Routs1  points 
for  0.600  inch  rods  closely  but  turns  sharply  up  at  a  volume 
ratio  of  unity  to  a  value  of  almost  38p  compared  to  Kouts0 
measurement  of  32. 5 • 

Since  Powell's  program  offers  a  choice  of  using  the  dif- 
fusion or  the  Fermi  formula  for  the  age,  for  our  purposes 
it  appeared  best  to  use  the  data  of  Kouts*  and  Chernick's  theory; 
since  Kouts0  data  was  obtained  from  the  physical  reality  by 
the  Fermi  formula,  the  machine  computation  will  also  be  based 
on  the  Fermi  formula. 

The  value  selected  was  a  constant  value  of  30  cm2  age  to 
indium  resonance  for  water  to  uranium  ratios  of  1  to  3 •   A 
1  cm2  correction  to  thermal  from  indium  resonance  gives  an 
age  of  31  cm2  for  all  enrichments  and  rod  sizes. 


Figure  F-I 

Fermi  Age  and  Migration  Area 
Measurements 
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Water  to  Metal  Ratio ,  V^/y 
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0.  Estimation  of  reflector  sayings 

Theoretically,  savings  for  an  infinite  water  reflector, ^ 
should  approach 

>"r 

Since  c^.  for  water  is  so  large,  it  will  overshadow  <^  for  g 

and  S3,  and 

Dc  2-  Dr 

At  room  temperature,  L  for  water  is 

T   =  —i ^  2.9 


/T 


However  experimentally,  values  range  from  approximately  six 
to  eight.   L  will  also  vary  with  temperature.   For  the  condi- 
tion in  this  thesis 


L  hot    r   /  ^   ^^    =  Nc   /  &Zl   <ZZT* 


x  1.35  ^  1-7 


1 
•  79 

We  assume  that  the  reflector  savings  will  vary  with  tempera- 


ture in  the  same  way  that  L  does 

Ah 


=  1.7 


By  selecting  the  end  of  the  experimental  range   closest  to  the 
theoretical  prediction 

Ah     =     1.7     Ac  =  1-7  x  6  =  10.2  =*  10  cm. 
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APPENDIX  0- 
OOMPUTATION  OF  CORE  RADIUS  AND  NUMBER  OF  RODS 

A  method  for  determining  neeessary  core  dimensions  for 
the  transfer  of  predetermined  thermal  powers,  given  the  unit 
cell  dimensions,  was  developed  from  Rohsenow,  Lewins,  and 
Barger' s  paper,  "Steady  State  Temperature  Distributions  in  a 
Nuclear  Reactor  with  End  and  Center  Fed  Coolant."  Rohsenow, 
Lewins,  and  Barger  develop  functions  P  and  Q,  for  the  tempera- 
tures of  the  fuel  centerline,  the  fuel  clad  interface,  the 
coolant  clad  interface,  and  the  coolant  bulk  temperature. 
Function  P  expresses,  in  dlmensionless  form,  the  temperature 
difference  between  the  point  in  question  and  the  ©oolant  in- 
let; function  Q,   expresses,  also  in  dlmensionless  form,  the 
heat  transfer  resistance  between  the  area  and  the  eoolant. 

It  is  shown  in  the  paper  that  the  maximum  P  along  a  chan- 
nel, and  hence  the  highest  temperature,  is 

p  2  =  Q2+l 

m 

Since  we  were  investigating  maximum  performance  of 
power  reactors,  the  limitations  placed  on  our  cores  were 
principally  those  of  temperature,  and,  of  these,  it  appeared 
that  temperature  at  the  fuel  rod  centerline  and  of  the  coolant 
bulk  temperature  would  be  the  most  critical.  Further,  it  ap- 
peared that  if  each  of  these  limitations  were  Just  met  at 
full  power,  that  the  most  efficient  reactor,  under  the  pre- 
scribed conditions,  had  been  designed. 

The  functions  given  by  Rohsenow,  Lewins,  and  Barger  for 

1.  Rohsenow,  Lewins,  and  Barger,  op  cit. 
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the  end  fed  case  are  as  follows: 
Coolant  Bulk  Temperature: 

P  -  7T^tc     x  wo     .    ,.ao     -  sin     JL    fo_ 
q«  A         af  2     af 

Q  ?  0 

t 

where  Atc     ~  t  coolant  max  -  t  coolant  inlet,  °F 

q"   =  heat  transfer  rate  at  midpoint  of  rod,  Ptu^ ig 

^m  *  »  hr  f  Xr 

w   =  weight  rate  of  coolant  flow,  lb./hr. 
c   a  coolant  specific  heat,  ftpoy 

A   =  channel  cross  section  area  f t 

a0  s  rod  heat  transfer  area,  midpoint  to  end  of  core,  ft2 

af  =  rod  heat  transfer  area  extended  through  reflector 
savings,  ftr 

Centerllne  Temperature: 


P  =  ir  Qt     x  wc  .  fo  _  gin  jr  ap_ 

q»      A      af        2   af 

=  ir  wc .  fo  ri+j.+R£X^  3 

*      A    af     h  Uc   R-o   ?fk  f 


E  (ft*©) 


where    h  =  film  heat  transfer  coefficient 

U0  s  effective  heat  transfer  coefficient  for  bond,  clad 
and  scale 

Ro  =  radius  of  fuel  region 

R  =  radius  to  outside  of  cladding 

1     -     Sfl  ln<Rc/Ro)     .lil 
ffc"~     "~^c"  *  hg  +Eb~ 

the   third  term  expresses  the  resistance  in  the  rod,   assuming 

that 

(f>    =     A  I0    (ft0R0) 

^tf-t  of  fuel  centerline  -  t  coolant  inlet 
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In  order  to  obtain  a  simultaneous  solution  for  given 

dtf8/]tc     determined  by  fuel  melting  and  coolant 
boiling  t's 


determined  by  materials  and  dimensions 
of  the  unit  cell 

in  each  case,  P2  was  not  equal  to  Q2-f  X.   In  the  equation  for 
the  coolant  bulk  temperature  thus  obtained ,  this  leads  to 
P  =  1,  and  the  resulting  equation  was  solved  for  w 

a, 


=   (1+ sin  5-^)  (A  q&) 


7T  4tc  .  c  .  ao/af 
Physically;,  this  means  that  we  have  set  coolant  flow  suffici- 
ent to  maintain  the  bulk  temperature  within  the  limits  set. 

This  value  of  W  was  then  substituted  in  the  fuel  center- 
line  equation;  that  is,  knowing  the  coolant  flow  rate,  we  de- 
termined the  amount  of  heat  that  could  be  generated  in  the  rod 
without  exceeding  the  temperature  limit  set. 

TT2    £*£    .    (1+Bln     ?    J£)    (Aq£)*    •    c2    .    (§£>2 

-2jrj*t    .    (1+slnf    gXAq")        o.foalnXao 
«*  T4t0.    o.    ao/af  A       «f  2     Sf 


+  sln2    T   !» 
2      af 


2  1 
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lm       -      l_  af  z  af  J 


2 


2 


where 


where   QB=    £  +  1    +  5*-  •     ^     '     !°£!2>JliT 

For  a  particular  series  of  reactors  In  which  the  rela- 
tion between  core  height  and  reflector  savings  remains  the 
same  and  in  which  the  same  temperature  limitations  apply, 
this  can  be  expressed  as 

«:=  T/ftR 

(l-Mln  5  ^J  2 

In  the  class  of  reactors  under  consideration,  core  heights 
will  be  of  the  order  of  8  feet  and  reflector  savings  of  the 
order  of  10  cm. 

Sln  ("?   a?)  =  8ln(J?  '  l|f)  =  sln  83"°  r  0.99^1 

a 

if  sin  2-  g°-  is  taken  as  constant  and  equal  to  1, 

T2  -     fat    (Atf     -  Jtc) 

q£  is  the  heat  transfer  rate  at  the  midpoint  of  the 
central  rod.   The  heat  transferred  in  the  entire  rod,  q0,  is 

q° s  q*  4°    cos  ( f  t? >  da 

*  a-       '    *  4     4af  qm 


=  qm  .   — £  sin  /if  ^o  )  -^ 
7T       l^  af' 


7T       v  *  af'      IT 
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And  the  average  power  developed  over  all  the  rods,  (assuming 
a  J0  flux  distribution)  q>  is 

q7TR2=  /q  J    2.405r.  27Tr  dr 

M  <   O   O    

7=-  s  a   „    J2 R»  j  /2./f05R  I 

where;      r    p  core  radius 

RO  a  R  +-/\ 

/4  =  reflector  savings 
The  thermal  power  generated  in  the  reactor, 

p  -  n  x  N      where  N  »  no.  of  rods 
*T     H 

JT  s„  =qox.83.|,  ^(^01*-) 

P   -   -83TTqo  R(R+A  )    t  /2.405R  \     A,,  =  cross  section  area 

T  ~ —   dl  1   R/ — )  u   of  unit  cell 

Au  v       ' 

This  may  be  solved  for  R,  and  the  number  of  rods  may  also 
be  found  from 

N  a   7TR2/AU 

Hot  channel1 factors  are  engineering  safety  factors  applied 
at  appropriate  points  in  the  heat  transfer  calculations  where 
uncertainties  exist  to  assure  a  safe  design.   They  may  also 
be  used  to  correct  for  the  non-uniformity  of  neutron  flux 
across  the  reactor,  but  this  was  accounted  for  in  the  process 
of  averaging  an  assumed  Bessel  function  J0  flux  distribution 
across  the  core  to  obtain  the  average  power  developed  in  one 


109. 


rod.   We  used  only  the  engineering  hot  channel  factors,  there- 
fore, and  adopted  a  set  very  similar  to  those  described  by 
J.  W.  Simpson,  et.  al.,  in  "Description  of  the  Pressurized 
Water  Reactor  Power  Plant  at  Shipping pe^rt,  Pa.M,  Geneva  Paper 
No.  815,  which  were  as  follows? 

?e        t      an  allowance  for  the  uncertainty  in  the  heat, 
transfer  coefficient,  h,  applied  to  reduce  h. 
prt   s   an  allowance  for  the  uncertainty  in  the  allowable 
heat  generation  rate  in  the  rods;  applied  to  re- 
duce q£  . 
F^^  s   an  allowance  to  assure  that  the  temperature 

rise  stipulated  is  sufficient  to  prevent  boil- 
ing, applied  to  reduce  A^q* 
The  values  of  these  factors  are  based  on  specified  tol- 
erances, known  ranges  of  uncertainty,  and  engineering  experi- 
ence and  judgment.   The  values  which  we  adopted  for  our  heat 
transfer  calculations  were  as  follows? 

*e     =  1-53 
Fq  a  1.08 

F2)t=  1-39 

By  analysis  similar  to  that  shown  above,  the  thermal 

fluxes  required  to  produce  burnout  and  local  boiling  can  be 

computed  and  compared  with  qJJ 

Burnout:    y     z-      "  F  ■<!  /  1  ,  l  1  2    - 
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.,.,,„,   qmLB,,         ^wo  (^*wo  °  FX? 

Boiling;   u,  ? —      s—t — r—^ a— ' 

*q  V  (l/h)^ 

Values  for  thermal  conductivities  and  the  film  coef« 

ficlent  were  selected  as  follows? 

Material  k     U/hr.ft. 


Uranium  oxide,  U02 

0.95 

Stainless  steel,  type  30^ 

10.6 

Zircalloy*-2 

8.15 

film  coefficient  h  =  6000  BtU/hr.ft.2  °F 

A  sample  calculation  of  core  radius  follows. 
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Sample  Calculations  for 
Core  Radius  U02  Fuel,  Stainless  Steel  Clad,  Core  No.  39 

T2   =   tf  {A*f  "   4tc)  =  *313  (^313-101)  =  18.2  x  10^ 
T    ^   ^270 

1     =    iiia  -  2.55  x  io-"-     f *  hr  F  -  2.55  x  10-* 

h        300O  Btu 

1         s       R°  ln  R°/R°     =      (-215)  1>    (:§&)  =  1.22  x  10-* 

V  kc  12  x  10.6 


fa 


o     =        . 211 


£c      x    _i__    x     Io(?fRo)-l 


E7 


f  TfRo  Io(^R0) 


..21' 


■  x       1     x   '0112       =  9I+.5  x  10~* 


-2^95^      72TT      T1061 


3ft 


98. 3  x  10"* 


q£        -     —       s  *2Z°  s  z^o.3  x  lO4 

m       "      Qr  I.OS  x  98.3  x  10~* 


Burnout   Che  ok 

1    4.  i.     =     2.55  x  10~*  +  1.22x10°*     -  3-77  x  10"* 

K  ^  u 

nM      -   2060   a  2060  =  547  xlO^  ^  2  q 

qmBo      £71      3.77x10-* 

Looal  Boiling  Check 


^m     :   ^0.3  x  10* 
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1/2  core  height,,   H  -  4  ft. 

Radius   to  outside  of     clad,        Rc     -      .215  in. 

Radius  of  unit  cell,   R-j     =   .880   cm. 

Co     -Qft*8ReH-     8(lH).3»10*j(.21?M»)  x  2,93  x  10-7 

~   .0676    mw 
R(H,10,xJlf|^)=     SUgg:  *  55X0 

Core  radius 9    R     =  9^.5  cm. 

R2 
Number  of  rods,  N  s  g~2  s  115|W 
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APPENDIX  H 
CORRELATION  OF  THE  POWELL  FOURIER  PROGRAM 
WITH  EXPERIMENTAL  DATA 

A  series  of  test  runs  was  made  with  the  Powell  Fourier 
program  for  Whirlwind,  using  data  from  experiments,  in  order 
to  compare  the  results  from  the  program  with  experimental 
results.   The  experimental  data  consisted  of  the  six  lattices 
tested  by  Krasik  and  Radkowsky— ^  (aluminum  clad  uranium  metal 
fuel  in  light  water,  with  characteristics  as  shown  in  Table 
H-I)  and  the  three  lattices  tested  by  J.  R.  Brown  et  al. 
(aluminum  clad  uranium  oxide,  UO2,  lattices  also  in  light 
water),  whose  characteristics  are  also  shown  in  Table  H-I. 

Since  values  were  published  for  the  fast  fission  factor 
and  the  resonance  escape  probability  for  these  lattices, 
these  values  were  used  in  the  computation  of  the  program 
input  M^ep10.   However,  values  were  also  calculated  for  g 
and  p  by  the  methods  shown  in  Appendix  D.   The  largest  error 
in  the  calculated  "^p"  was  1.9#,  and  the  average  about  \%. 
Adjusted  cross  sections  were  computed  as  discussed  in  Appen- 
dix C  and  in  the  Procedure.  At  first  these  were  adjusted 
to  an  arbitrary  density  and  atomic  weight  of  one,  and  this 
value  was  sent  into  the  computer  to  replace  the  values  for 
aluminum  which  were  a  part  of  the  program.   Later  it  was 
found  that  these  corrections  were  not  being  made  by  the 

1.  Krasik  and  Radkowsky,  op  cit. 

2.  Brown,  Nooroff,  Frantz,  Volpe,  and  Harris,  op.  cit. 
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machine,  and  the  cross  sections  were  then  adjusted  to  the 
aluminum  values ,  which  were  left  undisturbed. 

These  runs  gave  results  for  critical  mass  of  U-235  which 
were  consistently  too  high  "by  about  30# •   In  an  effort  to  find 
a  parameter  which  could  be  used  to  adjust  the  results  to  bring 
them  closer  to  the  experimental  results,  runs  were  made  in 
which  the  input  M'»£pM  was  varied  by  a  series  of  arbitrary 
factors.   This  is  the  input  in  which  the  greatest  uncertainty 
exists;  it  is  also  one  which  may  be  changed  without  requiring 
any  compensating  changes  in  other  inputs.   And  finally,  the 
calculations  are  quite  sensitive  to  it;   a  change  of  \%   in 
3?£  p  changes  the  critical  mass  by  about  2.5#»   The  results 
of  this  series  are  shown  in  Table  H-II  and  Figure  H-l  in  whioh 
fractional  variation  of  ygj>   from  the  calculated  value  is 
plotted  against  the  fractional  variation  of  the  resulting 
computed  critical  mass  from  the  experimental. 

It  will  be  seen  that  these  results  form  a  reasonable 
family  of  curves  and  appear  to  be  fairly  consistent.   In  order 
to  give  an  average  value  of  one  for  M25  (calculated)/!^ (experi- 
mental), an  arbitrary  factor  of  1.1  was  chosen  for  ^p/^ep 
(calculated).   This  method  of  correction  was  chosen  over  the 
alternate  method  of  using  lp£  p( calculated)  and  correcting  the 
resulting  M25  by  a  factor  of  3/1  »3  on  the  following  basis s 

1.  The  curves  were  steeper  at  7?£ p/7?£P  oalo.  -  1.1 
than  atftgp/9£p  calc.  s  1.0.   Thus  an  error  of 
1%   in  yev   would  make  an  M25  error  of  only  about 
2%s   against  an  error  of  5  or  6%   if  ip^p  were  in 
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the  other  range. 
2.  The  spread  of  l^calc./l^esp.  was  narrower, 
about  1  8*, at  *2?£p/^?£  p(calc.  )  =  1.1  than  that 
a^V£V/#£I>   calc.  s  1.0,  where  it  was  about 
I  10*. 
It  thus  appeared  that  both  the  error  and  the  range  of 
error  would  be  smaller  using  the  arbitrary  factor  of  1.1  in 
•»£  p  than  if  an  arbitrary  1/1.3  were  applied  to  critical 
masses  computed  from??£p  calc. 

There  is  some  measure  of  physical  justification  for  this 
correction.   In  relatively  tight  lattices  such  as  these, 
there  is  a  considerable  amount  of  epithermal  fission,  amount- 
ing to  as  much  as  10*  or  more  of  the  total  fissions,  acoord- 

1 
ing  to  Kouts.— -  If  these  did  occur  to  a  significant  degree, 

they  would  tend  to  reduce  the  mass  of  U-235  necessary  for 
criticality,  but  would  not  be  reflected  in  the  calculation. 
And  if  they  do  occur,  a  not  unreasonable  way  to  account. for 
them  seemed  to  be  by  a  fifth  factor  in  k^  ,  which  is  effec- 
tively the  correction  which  was  made. 


1.  Kouts,  Price,  Downes,  Sher,  and  Walsh,  op.  cit. 


116. 


I 


o 


^ 


117. 


TABLE  H-I 


Results 

of  Critical 

Experiments 

Core 

Vi/Vo 

fjf) 

Enr. 

Rod 
Dlam. 

Nr. 
Rods 

M(25) 
x*m°4 

P 

€ 

f 

La 

3 

1.3 

0.6M 

266 

1.44 

.850 

1,047 

,849 

B 

2 

1.3 

0.6W 

335 

1.82 

.794 

1.061 

,890 

C 

1.5 

1.3 

o.6M 

478 

2.58 

.720 

1,072 

,918 

D 

3 

1.15 

D.** 

380 

1.80 

.839 

1,050 

,836 

E 

2 

1.15 

0.6" 

463 

2.22 

.792 

1.065 

.881 

F 

3 

1.3 

.3871 

631 

1.42 

- 

- 

2. 
"12 

3 

1.3 

0.6M 

1268 

2.42 

.810 

1,039 

,875 

13 

4 

1.3 

0.6M 

1025 

1.98 

.835 

1.033 

.838 

14 

5 

1.3 

0.6" 

986 

1.90 

.870 

1.029 

.810 

1.  Cores  A,  B,  C,  D,  E,  and  F  are  uranium  fuel  measurements 
from  Krasik  and  Radkowsky. 

2.  Cores  12  „  13,  14  ar©  U0£  fuel  from  critical  measurements 
"by  Brown,  Noordof  f ,  Frantz,  Volpe  and  Harris  and  f ,  p,  g 
measurements  by  Krantz,  Smith;,  Klein,  and  Baer. 


TABLE  H-II 
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Comput 

ed  Results  for 

Varied  -Tp£p 

Core 

P 

M(25)oalc. 
(gm)   x  10^ 

A 

2.05 
1.95 
1.85 
1.75 
1.65 

1.11 
1.26 
l.kk 

1.69 

2.05 

B 

1.65 

2.82 

1.75 

2.28 

1.85 

1.91 

C 

1.73 

2.36 

1.83 

1.98 

1-93 

1.71 

D 

1.73 

2.68 

1.83 

2.2? 

E 

2.00 
1.875 

1.89 
2.6k 

1.755 

2.71 

1.625 

3.51 

1.50 

iK.88 

F 

1.71 

2.06 

1.81 

1.71 

1.91 

1J* 

12 

1.75 

3.4-0 

1.85 

2.93 

13 

1.95 
1.80 

2.57 
2.73 

1.90 
2.00 

2.37 
2.10 

1* 

1.86 

2.52 

1.96 

2.21 

2.06 

1.98 
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APPENDIX  J 

COMPUTATION  OF  THE  ADDITIONAL  ENRICHMENT 

REQUIRED  TO  ATTAIN  SPECIFIED  CORE  LIFE 

In  the  body  of  this  thesis,  a  series  of  cores  are  com- 
pared which  are  designed  to  have  equal  power  outputs;  that 
is,  they  are  capable  of  the  same  rate  of  heat  generation. 
In  order  to  select  a  true  optimum,  however,  the  designs 
considered  must  be  subject  to  a  second  requirement;   they 
must  be  capable  of  the  same  total  amount  of  heat  production, 
or  "power-time"  of  operation.   The  limitation  on  power-time 
is  that  the  reactor  must  remain  critical  throughout  the 
period  of  operation. 

The  principal  nuclear  reactions  affecting  the  reactivity 

of  a  low  enrichment  uranium  fueled  reactor  are  the  following? 

U-235  —  (n,f)  FP 
1 
1 (nfT)  u-236 

U-236  —    (n,    2fl  )  Pu-239  —    (n,f)  FP 

• (n,Y)  Pu-2^0  —    (n9Y)  Pu-24l  —    (n,f)  FP 

Production  of  the  fissionable  isotopes  Pu-239  and  Pu-24l  tends 

to  increase  the  reactivity,  and  their  destruction  to  decrease 

it  as  does  destruction  of  U-235.   Production  of  the  neutron 

absorbing  isotopes  U-236  and  Pu-240  and  the  neutron  absorbing 

fission  products  (FP  in  the  above  expressions)  tends  to  decrease 

reactivity.   Since  these  production  rates  will  be  different 

in  every  reactor,  depending  on  the  quantities  of  U-235  and 

U-238  present,  the  flux  at  which  the  reactor  is  operated,  and 
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the  resonance  absorption  probability  in  U-238,  each  Just 
critical  reactor  will  have  a  different,  and  relatively  short, 
allowed  time  of  operation  before  ultimate  loss  of  criticality. 
In  general,  it  will  be  necessary  to  increase  the  initial  en- 
richment to  a  value  greater  than  that  required  for  Just 
reaching  Initial  criticality  in  order  to  attain  practical 
power-time 8.   This  appendix  develops  a  method  for  finding  the 
factor  by  which  the  minimum  enrichment  for  initial  critical- 
ity of  a  given  reactor  must  be  increased  to  attain  specified 
power-time. 

While  the  core  life  is  expressed  in  power-time,  it  is 
useful  also  to  define  flux  time  as  a  measure  of  the  degree 
of  exposure  of  the  fuel.   The  number  of  reactions  which  take 
place  in  a  medium  of  given  ^  which  is  exposed  to  a  flux  f6 
for  a  time  T  is 

Since,  in  the  practical  case,  flux  will  be  varied  over  a 
period  of  time  due  to  changes  in  power  demand,  the  total  num- 
ber of  reactions  over  such  a  period  will  be 

T 


<£   0     0   (T)  dT 


The  quantity  /   T  0    (T)   dT  is  defined  as  the   flux  time,    t. 

Jo 

Following  the  method  of  Spinrad,  Carter,  and  Eggler, 
the  following  nomenclature  is  defined; 


1.  Spinrad,  B.  I.,  Carter,  J.  C,  and  Eggler,  C,  Reactivity 
Changes  and  Reactivity  Lifetimes  of  Fixed  Fuel  Elements 
in  Thermal  Reactors",  paper  no.  835  presented  at  the  Geneva 
Conference  on  the  Peaceful  Uses  of  Atomic  Energy,  1955* 
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£^28     =  N28  <3~28     eff 
^.25    =N25  (T25 

£  z+9     =  N49  ^  49 

^         =  N       (T 
£,   FP         FP       FP 


where  subscripts  25,  28, 
49,  and  FP  represent  U-235, 
U-238,  Pu-239,  and  pairs 
of  fission  products  respec- 
tively.  N,  the  number  of 
atoms  per  unit  volume  of 
the  reactor,  may  be  a.  func- 
tion of  flux  time. 


IOR,  the  Initial  conversion  ratio,  is  the  ratio  of  production 
of  new  fissionable  atoms  to  the  destruction  of  U-235  at 
startups 


ICR 


_.  w28  ^28  eff  =  SS£    +*1M1     ^    25  £   P^l-P) 


•25  u  25        ^.25    ^25 
where  P..  is  the  non-leakage  probability  during  slowing  down 
to  the  resonance  energy.   Since  the  effects  of  U-236,  Pu-240, 
and  Pu-24l  are  small  they  are  neglected  for  simplicity. 

For  the  low  enrichment  case,  differential  equations  may 
be  set  up  for  the  flux  time  changes  in  the  macroscopic  cross 
sections? 


£,. 


28 


o 
28 


jl£ 


dt 


25 


*  £ 


25        '  25 


Ji  4^9 
dt 

d  <h  FP 
dT 


cr 


-    cr 


49    *  ^'28 


<S  ?ft  "     2j 


49) 


FP 


£1 


1  +  <K 


■+ 


25 


^■49 


49 


Where  0\    is  the  ratio  of  capture   to  fission  cross   sections 
and  the   superscript   "o"   indicates  the  value   at   startup.      Solu- 
tions of  these  equations  are  obtained  as  follows? 


(1) 
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£25=  £.25  «-ff25t  <2> 

X,9  =     ^28   (X-."01***)  (3) 


^FP  = 

5  °    cr 

-4.25         FP 

i  +<     cr 

25      25 

Since 

^l.28 

ICR  = 

l-e 


1  +*49L        ^49 


-^ZS*     1      <^28   kr„t.   °FP   /,  49* 


l-e 


n(4) 


^Tz8  =     ICR  x    ^25  (5) 


And   equation    (k)  may  be  written  as 

illE=  _i £j?     i-e-^S*  +  ICR      fr„t  .   OJP  fi  e-°"*9* 

4  25       1+0<25     ^25  1  tc<4sL  ^9 


(6) 


or 


fe  =  -i-   1-e"  °"25*  +  ICR. 
M25  1  +^25  1   +  <X*9 


^9 


(7) 


Equation  (7)  expresses  the  total  number  of  fission  product 
pairs  produced  per  atom  of  U-235  initially  charged,  at  any 
flux  time,  t,  and  thus  the  total  number  of  fissions,  which 
is,  in  turn,  directly  proportional  to  the  total  heat  produced. 

The  productivity,  X  ,  is  defined  as  the  difference  in 
rates  of  neutron  production  and  loss  per  unit  of  flux  time. 
When  IT  =  o,  the  system  is  Just  critical. 

For  each  set  of  core  parameters  studied,  the  concentra- 
tions of  all  elements  represented  in  the  just  critical,  or 
"reference",  condition  may  be  found  from  the  cell  character- 
istics and  the  critical  mass  of  U-235  as  determined  by  the 
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machine  computation.   For  these  reference  systems  (denoted 
by  the  subscript  MRM),  productivity  is  zero. 
T  -  0  ■  &r  ■  CJ25-1)  ^25R  ~^28R  -^1R  "  LeakageR  =0  (8) 

After  the  required  flux  time  of  operation,  TT  must  again  be 
zero. 

&     -    (^5-1)  ^25  +  0]  49-D  <C*9  -  ^28  '  ^rL^ka8e- ^FP=0   (9) 

If  the  enrichment  of  the  reference  reactor  is  Increased 
by  a  factor  X 

W     =  (  25"1)  X  ^25r  "  ^28r  '  $1R  "  Leak«8eR  (10) 

Subtracting  equation  (8)  from  equation  (10) 

SV=  (^25-i)  ^25R  (X-l)  (ID 

Assuming  that  ^28   ^1  and  Leakage  do  not  change  appreci- 
ably during  the  period  of  operation,  equation  (10)  may  be 
subtracted  from  equation  (9): 

P|25-D(  %zrx  %25r)  -  0^9-1)  %k9  -  ^FP   -  -  Sir  (12) 

£\  FP  as   expressed  in  equation    (6)  may  be  broken  up   into  two 
parts 

^FP,  *     i£        ^JL     (l-e^25t)  (13) 


iK25     cr25 


.     IOR  ^25 


C 1, 0t 


rcrTOt      cr„    /,  .-  u*9 


^FP2   =     ^  ^      °  FP*  -      °_FP      (  1-, 
lt^49       L  % 

Substituting  values  from  equations    (2),    (3),    (5),    (13),    and 


(14) 
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(14)  into  equation  (12): 

en 25-D (x  ^25r  «"G'25t-x  ^5r).  lii^  ^r  (1.r^) 

1  R  1  tc<25     (r25 

4    (^49-D   X^25r     lCR(l-e-^9t) 


-  ICR 


C>Uo-l)<l-*e 


Rearranging  the  terms  and  dropping  out  the  common  factor  2|  2<R 

xTd-e-  ^Wl-W  _1     °"pp  U  I0R  f-  ^ 


)(1~W   I 

1    1+<^25  ^25^      (   1*K49 

=-X(^25-1)  +  (7  25-1)  (16) 

Equation  (16)  gives  the  required  factor  of  additional  enrich- 
ment as  a  function  of  flux  time  and  Initial  Conversion  Ratio. 

Appropriate  values  of  the  constants  in  equations  (7) 
and  (16)  for  the  conditions  specified  for  the  studies  are? 

FP  =  55  "barns 
^"25  =  431  0(   25  -  -204       77  25  -  2.0  8 

^9  =  687  <X  49  3  -^20       ^  49  =  2,°3 

Substituting  these  values,  numerical  forms  of  equations 
(7)  and  (16)  are  obtained: 
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x  ,  h™ . (17) 

1.08-1.19(l-e-43:Lt)t  IOR  [l.09(l-e-687T)-.0387t] 


^^—^l-e^431*)*  ICR[.30Ut-.44l(l-e--687*)J         (18) 


Nfp 


25 

Where  t  Is  flux  time  x  10"21 

Using  these  equations,  X  and  _Z^L  may  he  plotted  against 

N25 

flux  time  for  a  series  of   ICR's,    as   in  Figure  1. 

Since  hoth  ICR  and  Npp/N25  dep©nd  °n  X,  in  order  to 
determine  X  from  Figure  J-l  an  iterative  process  is  necessary. 
The  .steps  are; 

1.  For  the  specified  total  number  of  fissions  and 

from  the  critical  mass  of  U-235  for  the  reference 

/   o 
reactor,  compute  Npp/N25- 

2.  From  the  unit  cell  geometry,  enrichment,  fast 
fission  factor,  resonance  escape  probability, 
and  slowing  down  non-leakage,  compute  ICR  as 
in  equation  (1). 

3.  Enter  the  upper  curves  of  Figure  J-l  with  • 
Npp/N25R  and  *CRr,  and  find  the  flux  time  of 
operation  needed. 

4.  Using  this  flux  time  and  ICRr,  find  X  from  the 
lower  set  of  curves.   This  is  the  first  approxi- 
mation. 

5.  In  the  second  approximation,  enter  the  upper 
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set  of  curves  with 


NFP/N25       "  1/X(NFP/N23R) 


and       ICR 


hf^)      <28  \77l?<  £  Pid-P) 


These  will  give  a  different  flux  time  and,  from 

the  lower  curves,  a'  new  X. 
6.  This  process  is  continued  until  the  solution 

converges.   This  will  usually  require  three 

or  four  approximations. 
The  solution  can  be  somewhat  shortened  by  replotting 
the  data  as  in  Figure  J-2.   Here  flux  time  has  been  eliminated 
between  the  two  curves  in  Figure  J-l.   In  addition,  the  value 
of  Npp/N2i  plotted  for  each  value  of  X  has  been  modified  by 
the  value  of  X  so  that  Npp/N2eR  can  be  read  directly.   In 

finding  X  by  Figure  J-2,  Npp/^c  is  found  for  the  reference 
design  and  establishes  a  vertical  line  which  needs  no  fur- 
ther modification.   Along  this  line,  ICRr  can  be  used  to 
find  an  initial  estimate  of  X.   This  X  is  then  used  to  modi- 
fy ICR  as  in  step  (5)  above,  and  the  process  repeated  until 
convergence.   Figure  J-l  is  a  somewhat  clearer  representa- 
tion of  the  logical  steps  in  the  determination  of  X,  but 
Figure  J-2  will  give  slightly  quicker  convergence  than 
Figure  J-l,  eliminates  the  modification  process  on  Npp/N2?  > 
and  concentrates  the  data  in  one  set  of  curves. 
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